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Niobium  alloys  with  improved  resistance  to  oxidation  are 
needed  for  use  as  structural  alloys.  The  feasibility  of  forming 
protective  alumina  scales  in  Nb-Al-X  alloys  by  the  selective  oxidation 
of  aluminum  has  been  demonstrated  and  also  has  been  shown  that  a 
low  aluminum  diffusivity  is  the  major  factor  limiting  the  selective 
oxidation  of  aluminum.  Additions  of  certain  elements  such  as 
titanium  are  thought  to  have  a  beneficial  action  on  accelerating  this 
diffusivity  but  information  on  the  effect  of  alloying  additions  is 
limited. 

In  this  thesis,  diffusion  studies  in  the  binary  Nb-Al  system  by 
the  Matano-Boltzmann  method,  and  diffusion  analyses  in  the  ternary 
Nb-Al-Ti  system  through  concepts  such  as  Kirkendall  shift,  crossover 
displacements,  penetration  tendencies  and  composition  paths,  were 
conducted  to  evaluate  the  diffusiv'iy  of  A1  and  the  effect  of  Ti  on  Nb- 
based  alloys. 

Nb-Al,  Nb-Ti,  Nb-Al-Ti  (b.c.c.)  solid  solutions  were  prepared  by 
arc-melting  and  were  homogenized  in  a  high  vacuum  furnace  before 
the  fabrication  of  diffusion  couples.  Bonding  of  the  diffusion  couples 
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was  obtained  by  hot  pressing  in  a  temperature  range  of  700°  to  1000 
°C  for  two  hours.  A  fine  dispersion  of  yttiria  chopped  fibers  was  used 
as  markers.  The  diffusion  treatment  temperatures  were  varied  from 
1350°  to  1550  °C  and  1200  °C  to  1400  °C  for  the  binary  and  ternary 
studies  respectively.  Microprobe  analysis  was  employed  for  the 
evaluation  of  composition  profiles  in  all  the  diffusion  couples. 

Tne  interdiffusion  coefficient  for  the  ^  solid  solution  of  Nb-Al 
system  varied  linearly  with  A1  concentration  in  the  range  of  1-3  at% 
and  it  followed  a  linear  Arrhenius  behavior  as  a  function  of 
temperature  in  the  range  of  1350  °C-1550°C.  The  average  activation 
energy  for  diffusivity  of  A1  in  Nb  was  calculated  to  be  78.4  Kcal/mole 
with  a  frequency  factor  of  4.5  x  10*2  cm^/sec.  These  values  are 
smaller  than  those  reported  in  the  literature  for  a  higher 
temperature  range.  This  discrepancy  is  attributed  to  the  differences 
in  the  evaluation  method  and  the  temperature  range.  The  influence 
of  divacancy  mechanism,  which  is  known  to  operate  at  higher 
temperatures  in  Nb,  may  account  for  the  higher  diffusion  rate  and 
activation  energy. 

The  composition  path  and  penetration  tendencies  in  the  Nb-Ti- 
A1  system  suggest  that  Ti  is  the  fastest  element  in  the  P  solid 
solution.  Qualitatively  the  penetration  tendencies  correlate  with  the 
melting  point  of  the  alloys.  The  results  of  Kirkendall  shift 
measurement  were  found  to  be  in  agreement  with  the  calculated 
crossover  displacement.  The  results  obtained  in  this  study  suggest 
that  the  addition  of  Ti  to  the  binary  alloy  increases  the  general 
diffusivity  of  the  system. 
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CHAFrERl 

INTRODUCTION 


During  the  last  decade,  the  necessity  of  both  higher 
performance  requirements  and  new  applications  for  spacecrafts  has 
renewed  the  interest  in  the  development  of  advanced  aerospace 
materials.  The  aerospace  industry  has  set  the  reduction  of  aircraft 
weight,  extension  of  flying  ranges,  increment  of  payloads  and 
shortening  of  take-off  requirements  among  its  permanent  objectives. 
The  accomplishment  of  these  goals  demands  the  development  of 
materials  with  improved  physical  and  mechanical  properties  such  as 
strength  to  density  ratio,  fatigue  resistance,  high  temperature 
rupture  strength,  creep  resistance,  high  temperature  corrosion 
resistance  and  fracture  toughness.  Currently,  it  is  thought  that 
aircraft  engines  work  at  half  of  their  potential  horsepower  due  to 
limitations  such  as  strength  and  melting  point  of  the  existing  high- 
temperature  Ni-based  superalloys,  which  are  the  major  structural 
material  used  in  military  aircraft  turbine  engines.  State-of-the-art 
superalloys  are  limited  in  temperature  capability  to  1100  °C  . 
However,  there  is  an  immediate  need  for  a  metallic  material  capable 
of  withstanding  temperatures  up  to  1400  ®C. 

Because  of  their  special  properties  at  high  temperature, 
refractory  alloys  were  considered  in  the  development  of  aerospace 
materials  in  the  early  '60’s,  when  aircraft  and  aerospace  nuclear 
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propulsion  were  developed.  Among  these  refractory  metals,  the 
lower  density  niobium  alloys  were  investigated  for  turbine 
applications.  Their  long-term,  high-temperature  strength,  ductility, 
and  creep  resistance  are  well  retained  at  high  temperatures, 
comparable  to  those  obtained  in  nickel-base  superalloys,  but  their 
poor  oxidation  resistance  in  environments  rich  in  atmospheric  gas 
species  [1,  2]  and  the  arrival  of  superalloys  properly  cooled  halted 
their  development  [3].  At  present,  niobium  alloys  are  used  in  non- 
critical  sections  of  aircraft  engines,  like  exhaust  nozzle  components, 
where  their  high  melting  temperature  are  needed. 

A  recent  approach  to  solving  the  problem  of  poor  oxidation  of 
many  high  temperature  structural  materials  is  the  addition  of 
aluminum  as  an  alloying  element.  Alumina  AI2O3  is  a  stable  oxide 
that  usually  shows  a  good  adherence  to  the  metallic  alloys.  There  are 
two  main  problems  associated  with  the  oxidation  behavior  of  the 
aluminum  containing  high  temperatures  alloys.  In  the  case  of  alloys 
such  as  Fe-Cr-Al,  the  fast  diffusion  of  cations  through  the  alumina 
scale  leads  to  the  occurrence  of  the  oxidation  process  at  gas/oxide 
interface  and  consequently  voids  are  formed  at  oxide/metal 
interface,  which  results  in  oxide  spallation  and  poor  adherence  of  the 
oxide  to  the  alloy.  In  the  case  of  Nb  alloys  it  has  been  suggested  [4] 
that  the  slow  diffusivity  of  aluminum  prevents  the  continuous 
formation  of  alumina  and  transition  oxides  are  rapidly  formed  and 
hence  undesirable  high  oxidation  rates  are  achieved.  Oxidation 
studies  on  Nb-Al  alloys  [5]  have  indicated  that  a  minimum  aluminum 
content  of  50  at%  is  required  to  obtain  a  continuous  alumina  scale. 
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The  addition  of  alloying  elements  such  as  Ti,  Mo  and  V  has  been 
shown  to  reduce  the  oxidation  rate  of  Nb  alloys. 

Although  the  poor  oxidation  resistance  of  the  A1  containing  Nb 
alloys  has  been  related  to  the  low  diffusivity  of  Al,  there  have  been 
very  limited  studies  [6]  in  which  the  diffusivity  of  Al  in  Nb  has  been 
measured  indirectly.  The  purpose  of  this  research  was  to  evaluate 
the  diffusivity  of  aluminum  in  binary  Nb-Al  and  ternary  Nb-Al-Ti 
solid  solution  alloys.  The  experimental  approach  used  in  this  study  is 
based  on  the  diffusion  couple  technique  [7]  with  post-diffusion 
treatment  evaluation  of  composition  proHles  by  electron  microprobe 
analysis  (EPMA).  Optical  and  scanning  electron  microscopy  (SEM) 
methods  were  employed  for  the  evaluation  of  the  microstructure. 
The  alloys  were  prepared  by  arc-melting  of  high  purity  metals.  The 
diffusivity  measurements  were  conducted  in  the  temperature  range 
of  1350  ®C  to  1550  °C  for  the  binary  case  and  1200  ®C  to  1400  °C  for 
the  ternary  one. 

The  first  part  of  this  thesis  approaches  briefly  the  theory  of 
diffusion  from  both  a  phenomenological  and  an  atomistic  point  of 
view.  In  addition,  a  mathematical  analysis  of  the  theory  is  presented 
as  it  provides  the  basis  for  the  Fortran-77  source  code  [8]  used  for 
the  evaluation  of  the  coefficient  of  diffusion  of  Al  in  the  binary 
system.  The  basis  for  the  experimental  techniques  employed  through 
this  work,  as  well  as  a  review  of  the  previous  work  done  in  the 
binary  Nb-Al  system  and  the  ternary  Nb-Al-Ti  system  are  also 
presented  in  this  chapter.  Chapter  3  simply  describes  the  materials 
used  and  the  experimental  procedures  followed  throughout  this 
research.  The  results  for  the  binary  Nb-Al  system  and  for  the 
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ternary  Nb-Al-Ti  system  are  presented  and  discussed  in  Chapters  4 
and  5,  respectively.  Results  in  Chapter  4  include  a  discussion  about 
the  experimental  difficulties  found  during  the  preparation  of  the 
binary  alloy  Nb-Al  used  in  the  binary  couples.  The  discussion  of  the 
ternary  diffusion  results  in  Chapter  5  is  based  on  a  phenomenological 
approach  of  basic  concepts  such  as  composition  paths,  Kirkendall 
displacements  and  penetration  tendencies.  A  summary  of  results  for 
binary  and  ternary  diffusion  is  given  at  the  end  of  the  respective 
chapters,  with  the  final  conclusions  of  this  research  presented  in 
Chapter  6. 


CHAPTER  2 
LITERATURE  REVIEW 


This  chapter  is  divided  into  three  sections.  First  the  literature 
concerning  the  Nb-Ti-Al  system  is  reviewed.  The  subjects  prevalent 
to  this  section  include  phase  stability,  diffusivity  and  oxidation  in 
ternary  and  the  constitutive  binary  systems.  The  second  part  of  this 
review  concentrates  on  the  mathematical  formalism  and  mechanisms 
of  diffusion  in  binary  and  ternary  systems  consisting  of  single  phase 
crystalline  materials. 


2.1  Nb-Ti-Al  system 

2.1.1  Phase  diagram 

This  section  is  devoted  to  establish  the  limits  of  both 
tempe  ature  and  composition  for  the  b.c.c.  P  phase  in  the  Nb-Al-Ti 
ternary  phase  diagram. 

The  binary  diagrams  composing  this  ternary  system  are 
presented  in  Figures  2. 1-2.3.  The  Nb-Ti  system.  Figure  2.1,  features 
an  isomorphous-type,  complete  solubility  of  both  components  in  a 
b.c.c.  phase  at  temperatures  above  882  ®C  at  the  Ti  side  [9].  This 
phase  field  extends  to  room  temperature  as  Nb  content  increases. 
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Figure  2.1 


The  Ti-Nb  phase  diagram  [9]. 
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Figure  2.2 


The  Al-Nb  phase  diagram  [10] 
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The  other  phase  present  is  the  h.c.p.  a-Ti  phase,  which  is  stable  at 
high  Ti  contents  and  temperatures  below  882  °C. 

Nb-Al  and  Ti-Al  systems  exhibit  intermetallic  compound 
formation,  namely:  a2-Ti3Al  (DO19),  y-TiAl  (LIq),  and  Tj-TiAls  (DO22)  in 
the  Ti-Al  system,  and  6-Nb3Al  (A15),  a-Nb2Al  (DgB).  and  Ti-NbAl3 
(DO22)  in  the  Nb-Al  system,  in  addition  to  the  terminal  solid 
solutions.  The  solubility  of  A1  in  Nb  at  room  temperature  is  not  very 
well  established.  The  phase  diagram  [10]  suggests  a  solubility  limit  of 
7.5  at%.  Maximum  solubility  of  A1  in  b.c.c.  Nb  is  23  at%  which  is 
reached  at  1960  ®C.  The  P  phase  in  the  Ti-Al  [11]  phase  diagram  is 
stable  only  at  elevated  temperatures.  The  solubility  of  A1  in  jJ-Ti 
increases  with  temperature,  reaching  a  48.5  at%  maximum  solubility 
at  1475  “C. 

The  ternary  Nb-Al-Ti  system  has  been  the  subject  of  several 
investigations  during  the  last  decade  as  a  result  of  the  recent  activity 
on  development  of  alloys  for  high  temperature  applications.  It  is  a 
rather  complicated  system  featuring  several  transition  equilibria. 
Perepezko  et.  al.  [12]  and  Kaltenbach  et.  al.  [13]  recently  have 
performed  a  systematic  evaluation  of  this  system  including  mappings 
of  liquidus  projections  and  determination  of  isothermal  sections  of 
the  system  at  1200  ®C.  Both  studies  show  excellent  agreement  with 
each  other  and  studies  made  by  Benderzky  et  al.  [14]  and  Jewet  et. 
al.  [15].  Examples  of  liquidus  projections  and  the  isothermal  section 
at  1200  °C  are  given  in  Figures  2.4  and  2.5.  The  liquidus  surface 
features  several  four-phase  equilibria  (I-IV)  with  the  subscripts  1, 
2,...  indicating  the  temperature  sequence  (1  higher  than  2).  In 
accordance  with  [12],  the  binary  invariant  reactions  are  represented 
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Isothermal  section  at  1200  °C  in  the  Nb-Ti-Al 
ternary  system  [12]. 


as  pi  (peritectic)  and  ei  (eutectic).  There  is  only  one  saddle  point  in 
the  system  which  is  indicated  by  an  open  circle  in  Figure  2.4.  The 
arrows  inscribed  upon  the  lines  show  the  directions  of  the  slopes  for 
the  univariant  equilibria  (liquidus  valleys)  toward  lower 
temperatures. 

It  is  observed  that  the  liquidus  valleys  decrease  in 
temperature  from  the  high  temperature  peritectic  points  in  the  Nb- 
Al  side  swinging  clockwise  towards  the  Al  corner  of  the  triangle.  The 
P -phase  surface  exhibits  the  highest  melting  points  and  it  covers  a 
big  range  of  the  total  area  of  the  liquidus  projection  surface.  The 
large  solubility  of  Al  in  p  phase  decreases  with  temperature.  This  can 
be  observed  by  contrasting  the  exteiision  of  the  P -phase  field  of  the 
liquidus  projections  with  their  respective  isotherm  at  1200  ®C.  In 
accordance  with  these  projections,  the  solubility  limit  for  Al  in  the 
b.c.c.  solid  solution  depends  on  both  the  temperature  and  the  Nb  (or 
Ti)  content.  As  the  temperature  is  raised  and  the  Ti  content 
increased,  Al  becomes  more  soluble. 

2.1.2  Diffusivitv 


In  general,  'normal'  diffusion  in  metals  is  characterized  by  the 
following  properties 

1)  Do  and  Q  are  temperature  independent  over  a  wide  range  of 
temperatures, 

2)  Do  lies  in  a  range  of  0.05  to  5  cm^/sec, 

3)  Q  (cal/mole)  is  approximately  equal  to  34  Tm  (melting 
temperature,  K), 


4)  The  diffusion  value  for  substitutional  impurities  is  within  an 
order  of  magnitude  of  the  self-diffusion  rate  of  the  host  atoms. 
In  accordance  with  diffusivity  investigations,  there  are  a  number  of 
b.c.c.  metals  which  show  anomalous  diffusion  behaviors.  They  do  not 
present  "normal"  diffusion.  For  these  metals  the  Arrhenius  plot 
presents  a  curvature  that  is  more  pronounced  than  that  seen  in  f.c.c. 
metals,  for  which  is  a  result  of  divacancies  contributions  [17].  As 
mentioned  earlier,  the  Ti-Nb  system  exhibits  a  complete  solubility  in 
b.c.c.  phase  at  elevated  temperatures.  Ti,  among  others,  is  one  of 
these  b.c.c.  metals  with  anomalous  behavior  [18].  Murdock  et.  al.  [18] 
found  that  In  D  vs  1/T  plots  for  diffusion  of  TH^  in  Ti  were  fitted  best 
by  a  line  of  continuous  curvature  (Figure  2.6),  which  indicated  an 
increase  in  Q  and  Do  with  the  temperature.  They  proposed  the  model 
that  Kidson  [19]  used  for  the  diffusivity  of  Zr  (b.c.c.)  to  explain  their 
results.  In  accordance  with  this  model,  there  are  two  competing 
volume  diffusion  mechanisms  operating  in  the  two  temperature 
regions:  intrinsic  and  extrinsic.  He  [19]  attributed  the  low- 
temperature  'extrinsic'  region  to  the  presence  of  a  temperature- 
independent  concentration  of  vacancies  due  to  the  presence  of 
impurities,  like  oxygen  atoms,  which  cause  the  concentration  of  this 
type  of  vacancies  to  be  much  greater  than  the  concentration  of 
'intrinsic'  or  thermal  vacancies  at  low  temperature.  The  total  curve 
can  be  obtained  by  the  sum  of  2  exponential  terms  representing  each 
of  these  contribution.  The  resulting  equation  obtained  by  Murdock  et. 
al.  [20]  for  the  Ti  case  is: 

DTi44=  3.58x10-^  exp(-31,200/RT)  + 

1 .09  exp(-60,000/RT)  cm^/sec 


(2.1) 
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Figure  2.6 
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Temperature  dependence  of  the 
diffusion  of  Ti^  and  V^8  in  titanium 
[18]. 
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the  Arrhenius  plot  which  they  assumed  was  caused  by  a  divacancy 
effect  at  high  temperatures  and  was  fitted  by  the  expression: 

D=  (3.7±0.3)  exp(-4.54  eV/kT)  + 

(8±3)  X  10-3  exp(-3.62  eV/kT)  cm2/s  .  (2.2) 

Pontau  and  Lazarous  [24]  experimentally  evaluated  the  diffusivity  of 
Ti  and  Nb  in  b.c.c.  Ti-Nb  alloys  at  different  concentrations,  which 
were  chosen  to  avoid  the  a/p  transition  and  to  obtain  the  metastable 
ca  phase  [25].  This  co  phase  as  discussed  earlier  is  thought  to  enhance 
diffusion  in  other  b.c.c.  materials  like  Ti  and  Zr  [20,  26].  They  [24] 
determined  simultaneously  the  diffusion  coefficients  for  Ti  and  Nb 
radioactive  tracers  over  the  temperature  range  of  950®-1511  °C 
using  lathe-sectioning  techniques.  Their  results  are  presented  in 
Figures  2.8  and  2.9  and  Table  2.2.  They  found  no  evidence  of 
anomalous  behavior  in  the  diffusivity  behavior,  except  for  low  values 
of  Do  and  Q  compared  to  those  reported  by  [27].  Their  Arrhenius 
plots  were  moderately  curved,  but  due  to  lack  of  enough  points  an 
accurate  two  exponential  fit  to  their  data  could  not  be  performed. 
They  did  not  found  any  evidence  that  supported  the  enhancement  of 
diffusion  at  low  temperatures  by  the  presence  of  co -phase  nor  via 
dislocations,  and  concluded  that  an  operating  normal  relaxed 
vacancy-divacancy  mechanism  was  consistent  with  their  results. 

Newmann  and  Tolle  [28]  recently  analyzed  the  published  data 
in  self-diffusion  of  b.c.c.  metals  by  assuming  a  two-defect  model  and 
various  one-defect  models.  They  found  that  Nb  self-diffusion  data 
could  be  fitted  accurately  by  the  2-exponential  model  proposed  by 
Eizinger  et.  al.  [23].  They  pointed  out  that  oxygen  contents  up  to  800 
p.p.m.  have  no  influence  on  the  diffusion  coefficient.  In  accordance 
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Figure  2.8  Measured  diffusion  coefficients  of  Ti  (circles) 

and  Nb  (triangles)  in  different  Ti-Nb  alloys  as  a 
function  of  temperature  [24]. 
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Figure  2.9  Diffusion  coefficients  of  Ti  (circles)  and  Nb 

(triangles)  as  a  function  of  composition  for  3 
different  temperatures  [24]. 
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Table  2.2  Frequency  factors  Do  and  activation 
energies  Q  for  Ti  and  Nb  diffusion  in  Ti  and  Ti-Nb 
alloys  (From  Ref.  [24]) 


Compoaitioa 

Tracer 

Of  (cm* /sect 

0  (keal/mole) 

TI 

T! 

{4.s4:S:S)*io-« 

31 .30  a  0.43 

Nb 

31.03*0.43 

TI 

35.61*0.53 

Nb 

(1.79:5:1;)*  io‘* 

38.22*0.30 

Tl 

(3.i5:l:S)*io-* 

41.96*1.24 

Nb 

(i.ia:5:3)*io’* 

47.33*1.61 

Tl 

(2JSI  :l:?l)*io"' 

39.02  *  1.03 

Nb 

(2.98  :}:^)*  19'* 

61.73*1.45 
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with  their  analysis,  .c.  metals  can  be  divided  in  three  groups: 
Group  IVb  (Ti,  Zr,  Hf)  where  the  metals  show  an  abnormally  high 
diffusivity  at  low  temperatures.  Group  V  and  VI  (refractory  metals) 
which  show  a  more  or  less  normal  behavior,  and  Group  I  (alkaline 
metals)  which  fall  in  between  the  previous  two.  This  classification  is 
better  observed  in  the  normalized  plot  in  Figure  2.10.  They 
concluded  that  one-defect  mechanism  is  improbable  for  both 
refractory  and  alkaline  b.c.c.  metals,  but  they  could  not  discuss  either 
the  nature  of  the  diffusion  mechanisms  that  operate  at  high 
temperature. 

Hu  An  et.  al.  [29]  measured  the  interdiffusivity  of  the  Nb-Ti 
system  at  low  temperatures  (300  ®-350  ®C)  using  X-ray  diffraction 
techniques  in  compositionally  modulated  structures  with  an  initial 
content  of  50  at%  Ti  and  50  at%  Nb.  They  determined  a  value  lower 
than  10*20  cm2 /sec  for  the  interdiffusion  of  the  system,  which  is  in 
agreement  with  the  diffusivity  results  obtained  at  high  temperature 
by  Ugaste  and  Zaykin  [30]  using  tracer  technique.  Ugaste  and  Zaykin 
[30]  also  investigated  the  concentration  dependence  of  the 
interdiffusion  coefficient  in  this  system.  Figure  2.11  shows  a 
comparison  of  their  results  with  those  from  other  authors  which 
prove  to  be  completely  inconsistent.  Figure  2.12  also  gives  the 
concentration  dependence  of  the  activation  energy  Q  and  the  pre¬ 
exponential  factor  Dq.  Despite  the  inconsistency  in  concentration 
dependence  of  D  values,  the  Q  values  seem  to  be  consistent  with  the 
values  given  in  Table  2.2.  In  general,  D  increases  strongly  with  the  Ti 
concentration,  which  is  consistent  with  the  concentration  dependence 
of  the  solidus  temperature.  Consistently,  the  activation  energy 


Concentration  dependence  of  coefficient 
of  mutual  diffusion  in  Ti-Nb  system: 
1-1400®,  2-  1300®,  3-  1200®,  4-  1100®, 
5-  1000  ®C.  From  [30]. 
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Figure  2.12  Concentration  dependence  of  activation 

energy  and  pre-exponential  factor  of 
mutual  diffusion  in  the  Ti-Nb  system 
[301. 
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increases  with  the  Nb  content. 

The  only  information  available  on  the  diffusivity  of  A1  in  Nb  is 
the  back-calculated  diffusivity  from  the  evaporation  rates  of 
aluminum  from  dilute  Nb-Al  alloys  conducted  by  Nikolaev  and 
Bodrov  [6].  In  accordance  with  their  experiments,  the  aluminum 
diffusivity  follows  an  Arrhenius  behavior  that  can  be  expressed  as: 

Dai=  450+/-10  exp[-(102730  +/-  7000)/rt)  cm2/sec  . (2.3) 

They  evaluated  he  thermodynamic  and  kinetic  characteristics  of  the 
evaporation  of  aluminum  present  in  niobium  as  an  impurity  (10'2- 
10'3  at%)  by  atomic  absorption  spectroscopy  in  an  argon  atmosphere. 
Their  proposed  method  to  measure  the  diffusivity  of  impurities  in 
the  solid  phase  is  based  on  measuring  the  rate  of  decay  of  the  atomic 
absorption  by  plotting  the  kinetic  curves  for  the  evaporation  process. 
They  found  that  at  a  given  temperature  the  concentration  of 
aluminum  vapor  over  the  sample  stays  close  to  the  equilibrium  value 
for  a  short  period  of  time  (30  seconds),  followed  by  a  steep  fall  of  the 
vapor  pressure  curve  due  to  a  depletion  of  aluminum  from  the 
surface  of  the  sample.  After  measuring  the  absorption  of  aluminum 
vapors  from  pure  aluminum  and  alloys  containing  a  known  amount 
of  it,  they  calculated  the  A1  activity  coefficient  and  its  temperature 
dependence.  From  these  experiments  they  were  able  to  determine 
the  kinetics  of  the  evaporation  process,  that  is,  the  motion  of  the 
atoms  in  the  solid  phase.  During  the  first  part  of  the  process  the 
transport  of  the  aluminum  atoms  from  the  surface  of  the  sample  is 
achieved  by  evaporation.  The  subsequent  decay  of  the  absorption 
curve  denoted  a  rate  limitation  on  the  availability  of  aluminum 
atoms  caused  by  the  diffusion  of  the  aluminum  in  the  niobium.  They 
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implied  that  the  diffusion  step  of  the  evaporation  process  obeyed  the 
expression  131]: 

Ns(t)/NoS=  exp  [y2]  erfc  (y)  (2.4) 

where  Ns(t)  and  Nq^  are  the  aluminum  concentrations  at  the  surface 
of  the  sample  at  the  time  of  the  measurement  and  at  the  beginning 
of  the  experiment  respectively,  and 

y2=  p2tDs  (2.5) 

where  is  the  diffusion  coefficient  to  be  calculated  and  P  is  an 
experimental  constant.  In  accordance  with  their  research  this  method 
offers  speed  of  operation,  continuous  observation  of  the  diffusion 
process,  and  application  to  the  solid  state  diffusion  of  alloy 

components  at  the  impurity  concentration  level  (<10*2  at.%). 

However,  the  method  is  not  direct  and  is  based  on  many 

assumptions. 

Raissanen  et.  al.  [32]  investigated  the  diffusivity  of  A1  in 

poly  crystalline  ion -implanted  a-Ti  in  the  temperature  range  of  600 
®C  to  950  °C  using  ion-beam  techniques.  They  found  the  following 
expression  for  the  coefficient  of  diffusion: 

D  (cm2/sec)  =  7.4x10*2  exp[-1.62  (eV)/kT]  .  (2.6) 

Pokoev  et.  al.  [33]  also  investigated  the  Al  diffusivity  in  a-Ti  using  x- 
ray  diffraction  techniques  at  700-850  °C  in  a  vapor  deposited  Al 
layer  on  Ti  foil.  They  found  that  the  temperature  dependence  of  the 
diffusion  coefficient  is  given  by  the  following  equation: 

D  (cm2/sec)  =  9.7x10*^  exp  [27,500  (cal/mole)/RT]  (2.7) 
Chaplanov  et.  al.  [34]  conducted  annealings  on  Ti-Al  films  at  200-500 
®C  temperature  range.  Below  300  ®C  interdiffusion  occurred  with 
subsequent  formation  of  Ti-Al  solid  solutions.  However,  above  this 
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temperature  intermetallic  phases  such  as  TiAla,  TiAI  were  formed. 
The  formation  of  a  TiAls  has  been  reported  elsewhere  [35-37].  Zhao 
et.  al.  [35]  studied  the  reaction  kinetics  of  Ti  films  on  A1  substrate  at 
450-550  °C,  and  found  that  the  formation  of  a  TiAla  layer  always 
occurs  at  the  Ti/Al  interface.  Rao  and  Houska  [36]  studied  the  inverse 
case,  i.  e.  Al  film-Ti  substrate,  and  found  the  same  TiAls  layer 
formed.  However,  when  the  film  was  initially  heated  up  to  900  ®C, 
there  was  a  release  of  Al  into  a-Ti  in  the  first  stages,  but  longer 
annealing  times  eventually  led  to  the  decomposition  of  the  TiAla 
layer  with  a  substantial  release  of  Al  into  the  substrate.  According 
with  their  results,  the  variation  of  the  coefficient  of  diffusion  with 
composition  was  less  than  one  order  of  magnitude  when  the 
composition  changed  from  0  to  20  at%  (Figure  2.13)  [38,  39].  Because 
of  the  limited  solubility,  they  could  not  calculate  accurately  the 
coefficient  of  diffusion  for  the  TiAls  phase,  but  it  was  estimated  that 
it  tended  to  be  one  order  of  magnitude  smaller  than  the  diffusion  in 
the  a-Ti  (Al)  phase  at  900  ®C. 

Liu  and  Welch  [40]  made  a  survey  of  the  diffusion  data 
for  aluminum  in  both  a-  and  |3-Ti.  Tables  2.3  and  2.4  show  the  pre¬ 
exponential  and  activation  energy  of  aluminum  diffusion  in  a-  and  P- 
Ti,  respectively,  and  Figures  2.14  and  2.15  present  the  Arrhenius 
plots  coefficient  of  diffusion.  The  Al  diffusion  data  in  a-Ti  show  poor 
agreement  as  they  are  scattered.  As  to  the  diffusivity  data  for  P-Ti 
phase,  again  the  differences  in  the  reported  values  for  Al  are 
considerable,  and  no  clear  explanation  is  offered  about  it. 

The  diffusion  of  interstitial  elements  in  Nb-Ti  and  Ti-Al  alloys 
has  been  extensively  studied.  Vasileva  et.  al.  [41],  for  example. 
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Figure  2.13  Interdiffusion  coefficients  vs  at%  for  the 
diffusion  of  A1  in  a-Ti.  The  present  data 
obtained  at  900°  and  635  °C  are 
compared  with  results  published  by  1) 
[381,  and  2)  [391.  (From  Ref.  [361). 
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Figure  2.15  The  diffusivities  of  O,  Al,  and  V  in  P-Ti. 

The  numbers  in  the  plot  correspond  to 
the  series  numbers  in  Tables  2.4  [40]. 
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determined  the  activation  energies  for  N  diffusion  at  600-800  ®C  in 
Nb  (21.2  kcal/mole),  and  Nb-10  at%Ti  (26-26.2  kcal/mole).  Lauf  [42] 
studied  the  effect  of  additions  of  substitutional  solutes  to  Nb  in  the 
oxygen  diffusivity.  He  found  that  additions  up  to  5  %  of  Ti  resulted  in 
a  decrease  of  the  O  diffusivity  in  as  much  as  one  order  of  magnitude. 
It  was  suggested  that  this  was  caused  by  trapping  of  oxygen  by 
substitutional  solute  atoms,  with  a  substitutional/oxygen  binding 
energy  of  0.6-0.8  eV  for  the  Ti  case. 

Regarding  another  interstitial  element,  Cannelli  and  Cantelli 
[43]  have  conducted  extensive  research  on  the  presence  of  hydrogen 
in  Nb-Ti  alloys.  They  determined  the  diffusivity  of  this  interstitial 
element  by  means  of  the  Gorsky  effect  [44].  In  their  research  they 
found  that  substitutional  Ti  acted  as  a  trapping  center  for  H  and 
lowered  its  mobility,  as  shown  by  Figure  2.16. 

2.1.3  Oxidation 

Niobium  has  a  high  melting  point  and  a  good  strength  at  high 
temperatures  but  oxidizes  very  easily  at  temperatures  above  650  ®C. 
The  main  oxidation  product  is  Nb205  which  exhibits  paralinear 
kinetics,  i.  e.,  the  oxidation  rate  of  Nb  for  short  times  follows  a 
parabolic  behavior  and  becomes  linear  at  longer  exposure  times.  The 
reason  for  this  transition  is  due  to  the  cracking  of  the  Nb205  scale 
under  the  developed  compressive  growth  stresses.  An  extensive 
study  [45]  on  the  effect  of  doping  on  the  cracking  behavior  of  Nb205 
has  led  to  the  conclusion  that  Nb  alloys  capable  of  producing  a 
surface  oxide  other  than  Nb205  to  separate  the  alloy  from  the 
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Diffusion  coefficient  as  a  function  of 
reciprocal  temperature  of  unalloyed  Nb. 
and  apparent  diffusion  coefficient  for  2 
Nb-Ti  alloys  [43]. 


Figure  2.16 
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corrosive  environment  are  required.  The  most  effective  alloying 
additions,  that  is,  those  that  produced  the  lowest  oxidation  rates, 
were  found  to  be  aluminum,  chromium  and  iron. 

Svedberg  [5]  studied  the  effect  of  alloying  to  improve  the 
oxidation  resistance  of  Nb-based  alloys.  He  has  reported  that  the 
lowest  oxidation  rate  was  observed  for  NbAla,  which  forms  a 
protective  continuous  layer  of  alumina  AbOsby  the  selective 
oxidation  of  Al.  Any  other  Nb-Al  alloy  with  an  A1  content  less  than 
75  atomic  %  was  not  able  to  fonn  alumina  scales  at  1200  ®C.  Although 
the  reason  for  this  behavior  is  still  not  clear,  it  is  known  that  the 
formation  of  a  protective  layer  generally  requires  that  the  oxide  of 
the  added  element  to  be  more  stable  than  any  other  oxide  of  the 
major  component  of  the  alloy.  Since  alumina  is  very  stable  and 
diffusivity  of  oxygen  through  its  film  is  rather  slow,  aluminum  is 
preferred  as  an  alloying  element  to  improve  oxidation  resistance. 

Generally,  the  parabolic  rate  constant  of  a  material  is  a  common 
estimation  of  the  maximum  temperature  of  its  useful  resistance  to 
oxidation.  In  accordance  with  Wagner's  Diffusion-controlled  theory  of 
oxidation  [46],  when  the  cations  are  the  more  mobile  species  across 
the  oxide  scale  on  the  substrate  surface,  the  parabolic  rate  constant 
can  be  expressed  as: 

(m/A)  =  k"  t  .  (2.8) 

where  A  is  the  area  over  which  reaction  occurs.  In  this  case,  k"  is  also 
referred  to  as  a  'scaling  constant',  and  was  derived  considering 
measurements  of  the  mass  increase  of  the  specimen  m  as  the 
reaction  parameter  to  follow  the  oxidation  reaction.  Usually,  a  weight 
gain  of  approximately  20  mg/cm^  in  1000  hours  exposure  time. 
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corresponding  to  a  parabolic  rate  constant  value  of  1x10'^^  g^cm  '^/s, 
is  a  useful  limit  to  measure  oxidation  resistance.  In  accordance  with 
Figure  2.17,  alloys  forming  alumina  scales  in  NiAl  compound  will  be 
useful  up  to  1400  °C.  When  AI2O3  is  formed  in  Nb,  the  maximum 
useful  temperature  for  such  alloy  will  be  also  1400  °C,  and  not  1500 
®C  according  with  the  Nb205  melting  point  .  The  reason  for  this 
decrease  is  the  formation  of  a  eutectoid  between  mixed  oxides  of  A1 
and  Nb  which  would  lead  to  very  fast  rates  of  oxidation.  Thus,  Nb-Al 
alloys  may  be  used  at  temperatures  up  to  1400  ®C  for  extended 
periods  of  time  when  they  are  protected  by  alumina  scales.  The  main 
problem  of  applying  the  protective  layer  method  is  that  the  selective 
oxidation  of  aluminum  can  take  place  not  only  by  forming  a 
protective,  continuous  surface  film  at  slow  oxidation  rates,  but  also 
by  an  undesirable  rapid  internal  alumina  precipitation.  According  to 
Meier  and  Perkins  [4],  the  problem  that  complicates  the  transition 
from  external  to  internal  oxidation  is  the  formation  of  transient 
oxides  of  the  substrate  material  during  the  oxidation  process.  When 
transient  oxides  such  as  Nb205  form,  the  alumina  nuclei  formed  at 
the  alloy  surface  will  be  continuously  removed  into  the  growing 
transient  oxide,  avoiding  the  formation  of  such  protective  film  [47]. 
Therefore,  the  formation  of  transient  oxides  such  as  Nb205  must  be 
suppressed  so  that  an  alumina  film  could  be  developed.  Following  the 
analysis  by  Wagner  [48]  for  the  case  where  these  transients  oxides 
are  avoided,  they  estimated  that  the  transition  from  external  to 
internal  oxidation  in  Nb-Al  alloys  occurs  at  a  minimum  content  of  A1 
in  the  alloy,  which  is  function  of  factors  such  as  diffusivity  of  A1  and 
O,  alloying  additions,  and  oxidation  temperature,  as  established  by 
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the  following  equation: 

N‘>ai  >  [  (^tg*/3)  •  No(s)  •  (DoV^)/(DaiVox)  ]1/2  .  (2.9) 

which  encourage  the  formation  of  the  alumina  at  lower  aluminum 
compositions  in  the  alloy  by  an  increase  in  the  aluminum  diffusivity 
and  a  decrease  in  the  solubility  and  diffusivity  of  oxygen.  Estimations 
made  by  Perkins  and  Meier  [4]  indicated  that  it  is  very  difficult  to 
form  a  protective  alumina  film  on  binary  Nb-Al  alloys  at  low  enough 
aluminum  content  to  consider  them  still  Nb-base.  In  accordance  with 
observations  from  other  authors  [49],  the  interdiffusion  coefficient  of 
b.c.c.  Nb-Al  alloys  has  been  found  to  increase  considerably  with 
aluminum  content,  and  the  value  reported  by  Bodrov  and  Nikolaev 
(section  2.1.2)  is  probably  an  underestimate.  It  has  been  suggested 
[SO]  that  the  probability  of  formation  of  an  alumina  external  layer 
would  be  enhanced  by  accelerating  the  diffusivity  of  aluminum  by 
means  of  the  use  of  other  alloying  elements.  In  accordance  with  some 
authors  [4],  an  increase  in  the  diffusivity  of  aluminum  of  about  2-3 
orders  of  magnitude  is  possible  if  an  extended  b.c.c.  solid  solution 
existed.  The  elements  that  increase  the  solubility  of  aluminum  in 
b.c.c.  Nb  are  titanium,  iron,  chromium  and  vanadium. 

The  solubility  of  oxygen  in  the  Nb-based  alloys  can  be 
suppressed  by  the  additions  of  alloying  elements  by  the  gettering 
effect  of  these  elements.  As  mentioned  in  the  previous  section,  Ti  is 
one  of  these  elements  which,  in  addition,  forms  an  oxide  that  has  a 
desired  stability  between  AI2O3  and  Nb  oxides.  As  explained  by  Lauf 
et.  al.  [51],  the  addition  of  1-5%  of  solutes  such  as  Ta,  Ti,  V,  or  Zr  to 
Nb-based  alloys  decreases  the  O  diffusivity  significantly. 
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The  topics  of  internal  and  external  oxidation  of  alloy  solutes 
have  been  studied  extensively  [50,  52,  53].  It  is  clear  from  these 
investigations  that  oxidation  protection  can  be  obtained  by  means  of 
complex  mechanisms  that  involve  the  diffusivity  of  O  into  the  metal 
substrate  to  form  oxide  particles  or  films  that  require  a  solute  to 
diffuse  out  such  substrate.  Titanium  would  appear  to  be  a  probably 
effective  addition  for  increasing  Dai  and  decreasing  Do  and  No  in  Nb- 
based  alloys.  Additions  of  20%  Ti  should  nearly  double  the  solubility 
of  aluminum  in  niobium  at  1400  ®C  as  it  forms  a  continuous  b.c.c. 
solid  solution  with  niobium  which  can  dissolve  42-48  at.%  of 
aluminum  at  1400  ®C  [4]. 

2.2  Diffusion  Theory  and  Mechanism 
2.2.1  Binary  Diffusion 

When  two  different  components  are  placed  in  contact  and 
allowed  to  interact  at  a  moderate  temperature,  thermodynamics 
predicts  that  a  variety  of  compositions  and  phases  may  form.  Under 
adequate  experimental  circumstances  the  interaction  between  the 
two  components,  so-called  Diffusion,  eventually  would  lead  to  the 
reproduction  of  all  the  phases  between  the  end  compositions. 
Diffusion  is  defined  as  the  movement  of  atoms  in  a  solution  as  a 
consequence  of  an  initial  chemical  potential  gradient,  usually 
expressed  as  a  composition  gradient,  which  the  thermodynamically 
unstable  system  try  to  eliminate  by  going  to  its  lower  energy  state 
(equilibrium).  When  the  interaction  between  the  two  components  is 
conducted  in  single  phase,  the  matter  will  flow  in  a  manner  that  will 
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.decrease  the  concentration  gradient  with  the  net  flow  of  matter 
ceasing  when  no  chemical  potential  gradient  exists  anymore  ,  that  is, 
when  homogeneity  in  the  system  is  reached.  Specifically  in  solids,  the 
diffusion  is  conducted  by  either  a  vacancy  mechanism  (substitutional 
solid  solutions)  or  interstitial  displacement  of  solute  atoms 
(interstitial  solid  solution)  depending  on  the  size  difference  between 
the  solute  and  solvent  atoms  making  up  the  solution.  Next,  a  brief 
analysis  of  the  fundamentals  of  the  theory  of  diffusion  for  the 
diffusivity,  mostly  in  substitutional  solid  solutions,  is  presented. 

2.2. 1.1  Pick's  Laws  of  Diffusion 

Adolph  Pick  [54]  proposed  in  1885  an  equation  that  fitted  the 
fact  that  the  flow  of  matter  of  component  i,  Ji,  across  a  given  plane 
goes  to  zero  as  the  system  reaches  equilibrium,  i.e.,  the  concentration 
gradient  vanishes: 

Ji  =  -  Di  [9ci/ax]  (i=  1,2)  .  (2.10) 

where  Ci  is  the  concentration  c  f  component  i,  and  x  is  the  distance 
taken  parallel  to  the  concentration  gradient  of  such  component.  The 
proportional  factor  in  equation  2.10  is  known  as  the  coefficient  of 
diffusion  (or  simply  diffusivity).  Originally  D  was  conceived  as  a 
constant  at  a  fixed  temperature,  which  later  on  was  proved  to  be 
true  only  for  gas  solutions. 

Equation  2.10  is  valid  only  when  the  considered  system  is 
under  steady  state  conditions  (no  concentration  changes  in  a  point 
with  time).  To  take  into  account  the  variation  of  composition  with 
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time.  Pick's  first  law  is  inserted  into  an  equation  describing  a  local 
conservation  of  mass.  'Iik.  Pick’s  second  law  can  be  expressed  as: 

dci/dt  =  divJi=  -div  [  -D  Vckl  .  (2.11) 

This  equation  has  the  properties  of  a  continuity  equation  and  stems 
only  from  the  conservation  of  matter.  There  are  two  classical 
solutions  to  Pick's  second  law,  which  are  discussed  here, 
a)  D  constant,  steady-  and  non-steadv  states 
If  D  is  independent  of  position  equation  2.11  can  be  expressed 


as 

dci/dt  =  D  V[Vck]  =  =  D  V2ck  .  (2.12) 

where  V2  is  the  Laplacian 

^2  =  02 /3x,3y,9z  . (2.13) 


if  a  steady  state  exists,  dddt  equals  zero.  Assuming  the  simplest  case 
(unidirectional  diffusion)  equation  2.12  reduces  to  D  0.  This 

case  does  not  occur  very  often.  When  the  diffusivity  is  not  a  function 
of  composition  (or  position  for  that  matter)  but  non-steady  state  is 
acting,  equation  2.11  is  expressed  as: 

dc/dt  =  D  02c  /9x^  .  (2.14) 

A  solution  to  this  equation  will  express  the  concentration  of  a 
component  i  as  a  function  of  position  and  time,  c(x,t),  under  the 
adequate  initial  and  boundary  conditions.  In  the  case  of  semi-infinite 
solids: 

c(x,t)  =  c’/2  [1  +  erf  (x/2VDt)  .  (2.15) 

where  c’=initial  solute  concentration  of  the  richer  side  in  component  / 

z 

erf=  error  function=  (2/^^)  J  (’U^) 

0 


The  coefficient  of  diffusion  varies  with  composition  in  real 
experiments.  The  dependence  of  D  on  composition  forces  the  Pick's 
second  law  to  be  written  as: 

53ci/5at=  -  3  [Di(aci/ax)]  /ax  .  (2.1 6) 

This  is  an  inhomogeneous  differential  equation,  which  its  solution  is 
extremely  difficult  in  the  best  case.  The  most  common  solution 
employed  in  metallurgical  studies  on  diffusion  is  the  so-called 
Boltzmann-Matano  analysis,  which  allows  D  to  be  calculated  from  a 
concentration  profile  plot.  This  method  is  based  on  a  solution  of  Pick's 
second  law  proposed  by  Boltzmann  [55]  in  1884  and  was  proposed 
for  the  first  time  by  Matano  [56]  in  1933.  After  a  diffusion  couple 
has  been  prepared  and  annealed  at  a  fixed  temperature  during  a 
time  t  the  experimental  concentration-penetration  curve  c(x,t)  is 
experimentally  determined  (e.g.  microprobe  analysis).  The  resultant 
curve  should  not  be  necessarily  symmetric  and  in  that  case  an  error 
function  solution  can  not  be  applied.  A  solution  to  equation  2.16  may 
be  obtained  by  applying  the  so-called  lamda  substitution  \  [57]: 

A.  =  x/Vt  .  (2.17) 

which  can  convert  the  partial  differential  equation  2.16  into  a  total 
differential  equation  relating  c  and  the  variable  X.  This 
transformation  results  in: 

[-X/2]  dc  =  d[D  dc/dA,]  .  (2.18) 

Equation  2.18  is  constituted  only  by  total  derivatives  since  c  is  only 
function  of  X.  Integrating  2.18: 
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c  c  c 

-I  X,/2  dc  =  |d  [D  dc/dx)  =  [D  dc/dXl  I  .  (219) 

c  -  c  -  c  - 

where  c-=original  composition  (constant) 

c=  variable  composition  at  point  A  (or  j 


=  D(c)  dc/dX  -  D(c-)  dc/dX  .  (2.20) 

Arranging  2.19:  D(c)  =  [  -J  X/2  dc]  [  dX/dcl  .  (2.21) 

Differentiating  equation  2.17  at  a  fixed  time  gives: 

dX  =  dx/t-l/2  .  (2.22) 

Substituting  in  2.21  yields: 

c 

D(c)  =  -l/2t  (dx/dc)c  I  X  dc  .  (2.23) 

c- 


This  solution  requires  a  proper  evaluation  of  the  x-0  point 

corresponding  to  the  original  interface,  which  may  be  evaluated  by 

applying  a  conservation  of  mass  condition.  A  procedure  to  locate  the 

original  interface  from  the  composition  profile  considers  Figure  2.18. 

The  total  area  under  the  curve  is  given  by: 
c 

j  c(x)  dx  .  (2.24) 

c  - 

The  plane  at  which  x  =  0,  so-called  Matano  interface,  is  determined  by 

finding  the  value  of  x  in  the  diffusion  zone  {xq)  such  that 
c  c  c 

J  c-dx  =  J  c+  dx  =  Jc(x)  dx  .  (2.25) 

c-  c-  c- 

Graphically,  the  Matano  interface  is  the  line  that  makes  the  two 
hatched  areas  in  Figure  2.18  equal.  This  figure  also  presents  a  sketch 
of  the  area  under  the  curve  {fxdc)  and  the  inverse  of  the  slope  dx/dc 
as  a  function  of  composition  (position).  Note  that  the  Matano 
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Figure  2.18  Typical  composition  profile  of  an  annealed 

diffusion  couple.  The  Matano  interface  is 
positioned  to  make  the  hatched  areas  on  (cross 
hatched  area),  and  the  tangent  values  employed 
to  calculate  the  intertdiffusion  coefficient  at  a 
composition  c=  O.I  cq. 
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interface  is  characterized  by  a  maximum  at  the  area  under  the  curve 
and  a  minimum  in  the  inverse  of  the  slope.  The  error  in  the 
calculated  D  is  largest  at  the  ends  of  the  curve  where  c-/c  + 
approaches  zero  or  one  as  in  those  regions  the  integral  is  small  and 
dxldc  large.  As  mentioned  before,  the  Matano-Boltzmann  method  is 
widely  employed  in  metallurgical  diffusion  studies  as  it  is  very 
useful  for  determining  D  over  a  range  of  compositions  in  alloys. 
Among  its  advantages  are  the  wide  availability  and  development  of 
both  electron  analysis  techniques  and  computers,  which  make  this 
procedure  highly  automated.  However,  one  detriment  in  the 
procedure  is  that  equation  2.23  assumes  that  the  atomic  volume  of 
the  alloy  is  independent  of  composition,  which  is  a  seldom  case. 

2.2. 1.4  Binary  Diffusion  bv  a  Vacancy  Mechanism.  The  Kirkendall 

Effect  and  Darken  Analysis 

Two  mechanisms  have  been  proposed  to  explain  diffusion  in 
substitutional  solid  solutions.  The  first  is  the  cooperative  movement 
of  atoms.  Zener's  ring  or  direct  interchange,  described  later  in  this 
section,  are  examples  of  this  mechanism.  However,  energetic 
considerations  and  the  fact  that  the  Kirkendall  effect  occurs,  ruled 
out  these  mechanisms  as  possible  ways  for  substitutional  diffusion  to 
occur.  The  second  kind  of  mechanism  implies  the  movement  of  a 
single  atom  at  a  time,  and  the  example  is  the  universally  accepted 
mechanism  of  vacancy  diffusion.  In  this  case,  the  atom  of  an  adjacent 
site  simply  jumps  into  the  vacancy,  but  for  that  to  happen,  a  small 
displacement  of  the  neighboring  atoms  is  required,  which  constitutes 
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one  of  the  energy  terms  required  for  diffusion  encountered  in  the 
expression; 

D=  a  a2  Z  u  exp*(^Gni  +  AGf)/RT  .  (2.26) 

where  Z=  coordination  number 
a=  lattice  parameter 
a=  crystallographic  geometric  factor 

o=  vibration  frequency  of  a  solute  atom  in  an  interstitial  site 
AGin=  free  energy  per  mole  needed  for  solute  atoms  to  jump 
between  interstitial  sites. 

AGf=  free  energy  term  per  mole  associated  with  the  formation 
of  vacancies. 

The  pre-exponential  factor  is  fairly  constant  with  the  temperature 
and  can  be  joined  in  a  single  term  called  frequency  factor  Dq. 

Knowledge  of  the  operating  mechanisms  of  diffusion  in  b.c.c. 
materials  is  not  as  extended  as  that  for  f.c.c.  metals.  Self-diffusion 
and  impurity  diffusion  in  b.c.c.  metals  occurs  via  vacancies.  In  a  b.c.c. 
lattice  the  barrier  for  the  jump  of  one  atom  into  a  neighbor  vacant 
site  is  more  complex  than  it  is  in  f.c.c.  structures  because  of  the 
saddle  point  barrier  with  double  maxima  in  energy  that  the  hopping 
atom  must  overcome  [58].  Fritzsch  et.  al.  [59]  studied  the  diffusion 
via  vacancies  as  the  dominant  mechanism  in  b.c.c.  metals  by 
calculating  the  parameters  of  a  given  interaction,  such  as  the 
influence  of  lattice  constant  and  binding  energy,  between  the  atoms 
and  the  lattice.  They  considered  three  types  of  jumps  as  exemplified 
in  Figure  2.19,  and  calculated  their  potential  curves  (Figure  2.20), 
which  show  that  the  jump  to  the  next  nearest  neighbor  type  A  is 
favored  due  to  the  smaller  distortion  in  the  hopping  event.  They 
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Figure  2.19  Atomic  jumps  in  a  b.c.c.  crystal 

considered  in  the  model  calculation 
made  by  [59].  a)  Type  A:  jumps  to  the 
nearest-neighbour  position,  b)  Type  B: 
jumps  to  the  second  nearest-neighbor 
position,  c)  Type  C:  jumps  to  the  third 
nearest-neighbour  position. 
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Figure  2.20  Potential  barriers  for  the  cited  types  of 

individual  jumps  starting  at  r=  0  and 
ending  at  r=  1.  The  maximum  value  of 
each  curve  corresponds  to  the  migration 
energy  for  the  particular  hopping  event 
[Z9]. 
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found  a  strong  proportionality  between  migration  energy  and  the 
binding  energy  of  the  crystal  for  the  self-diffusion  case  (Figure  2.21), 
which  is  in  good  agreement  with  published  experimental  data.  As  to 
the  impurity  diffusion  case,  it  was  concluded  that  the  bigger  the 
impurity  atom,  the  higher  would  be  the  migration  energy,  which  also 
agrees  with  experimental  data  presented  by  Tiwari  et.  al.  [60].  The 
binding  energy  found  by  [59]  has  only  a  small  effect  in  this  migration 
energy. 

Regarding  with  the  solution  to  the  diffusion  equations 
presented  before,  it  should  be  pointed  out  that  the  Matano- 
Boltzmann  solution  gives  only  one  diffusion  coefficient,  D ,  that 
completely  describes  the  process.  In  this  kind  of  experiment,  the 
coordinate  system  employed  is  fixed  relative  to  the  ends  of  the 
specimen,  and  it  is  observed  that  the  diffusivity  gives  rise  to  the 
movement  of  one  part  of  the  diffusion  couple  relative  to  the  other. 
Therefore,  ahy  movement  of  lattice  planes  relative  to  the  ends  of  the 
couple  is  recorded  as  a  flux  and  affects  D  (equation  2.10),  even 
though  such  translation  does  not  correspond  to  any  jumping  of  atoms 
from  one  site  to  another.  Since  the  process  takes  place  in  a  chemical 
concentration  gradient  the  resultant  D  represents  the  cnemical 
diffusion  coefficient,  also  called  the  interdiffusion  coefficient.  For  the 
complete  description  of  the  diffusion  behavior  of  binary 
substitutional  solid  solutions,  however,  at  least  five  different 
diffusion  coefficients  are  needed  for  each  concentration,  namely: 

D=  interdiffusion  coefficient 

Da,B=  intrinsic  diffusion  coefficients  of  A  and  B  components 

D*a,b=  tracer  diffusion  coefHcients  of  A  and  B  components 


Migration  energy  for  self-diffusion  in 
dependence  on  the  binding  energy  of 
the  crystal  in  comparison  with  the 
experimental  data  refered  therein.  The 
line  drawn  presents  the  calculated 
relation  by  [59]. 


The  tracer  diffusion  coefficients  arises  from  the  fact  that  studies 
showed  that  the  self  diffusion  coefficients  [61]  are  in  general 
different  for  the  two  elements.  In  self-diffusivity  studies  one  studies 
the  diffusion  of  a  radioactive  isotope  solute  in  a  matrix  of  the  same 
(non-radioactive)  metal.  Except  for  a  small  mass  difference,  the  two 
atomic  species  are  identical.  This  mass  difference  permits  to  calculate 
the  jump  freq’iency  of  the  isotope  and  a  self-diffusion  coefficient, 
which  for  b.c.c.  metals  is  given  by: 

Db.c.c.  =  a2/8t  .  (2.27) 

where  x  is  the  dwell  time  of  an  atom  in  a  lattice  position.  Thus,  a 
relationship  between  the  chemical  diffusion  coefficient  and  the  self¬ 
diffusion  coefficients  must  be  found.  Such  relation  was  brought  about 
in  1947,  when  Kirkendall  and  Smigelskas  [7]  conducted  diffusion 
experiments  on  diffusion  couples  formed  by  a  rectangular  bar  of  70- 
30  brass  plated  with  0.1  inch  of  copper.  The  original  feature  of  their 
work  was  the  incorporation  of  Mo  wires  in  the  interface  of  the  two 
materials  as  markers  (Mo  is  insoluble  in  either  Cu  or  Zn).  After 
annealing,  a  thin  slice  of  the  couple  was  cut  and  metallographically 
prepared  to  measure  the  distance  between  markers.  The  experiment 
was  repeated  at  different  annealing  times.  In  all  cases  was  found 
that  the  markers  distance  decreased  an  amount  proportional  to  the 
square  root  of  the  annealing  time,  even  after  subtracting  the  small 
volume  change  in  adding  Cu  to  brass.  The  only  way  to  explain  this 
markers  shift  was  that  each  of  the  two  atomic  forms  moved  with 
different  velocities,  i.  e.,  that  the  flux  of  Zn  atoms  outwards  across  the 
markers  plane  be  much  larger  than  the  flux  of  Cu  atoms  inwards 


across  the  same  plane.  Later  work  in  different  systems  prove  the 
generality  of  this  concept. 

The  Kirkendall  effect  can  be  taken  as  a  confirmation  of  the 
vacancy  mechanism  of  diffusion  in  cubic  materials  (f.c.c.  and  b.c.c.). 
As  mentioned  earlier,  another  possible  diffusion  mechanism  in 
substitutional  solid  solutions  of  b.c.c.  metals  is  the  ring  mechanism 
proposed  by  Zener  [62].  Here,  the  energy  needed  for  diffusion  is 
smaller  due  to  a  more  open  structure;  however,  experiments  on 
many  diffusion  couples  composed  of  a  number  of  different  b.c.c. 
metals  showed  that  a  Kirkendall  effect  occurred.  Interchange 
mechanisms  such  as  the  ring  one  require  that  the  rate  at  which 
atoms  of  A  move  must  equal  the  rate  at  which  B  atoms  move  in  the 
other  direction.  Consequently,  no  Kirkendall  effect  would  be 
observed. 

Kirkendall's  experiment  showed  that  the  atoms  of  two  species 
in  a  diffusion  couple  move  at  different  rates.  Thus  in  accordance  with 
this  concept,  the  interdiffusion  coefficient  previously  defined  should 
be  used  only  to  describe  the  flow  of  components  relative  to  some 
component  averaged  convective  velocity  and  basically  is  the  same 
for  components  A  and  B.  Kirkendall  experiment  showed  that  using 
the  right  reference  frame  for  the  motion  of  the  atoms,  one  can  get  a 
different  result  for  the  diffusion  behavior  of  each  of  the  components 
and  talk  about  a  behavior  which  is  intrinsic  to  each  of  the 
components. 

In  1948,  Darken  [63]  published  an  excellent  phenomenological 
analysis  of  diffusion  of  alloys  based  on  such  experiment.  He 
contended  that  the  markers  motion  was  due  to  an  imbalance  of  flows 
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of  the  two  components,  and  that  each  component  had  a  flow  of  rate 
of  his  own,  bringing  along  the  definition  of  intrinsic  diffusion 
coefficients.  Consider  Figure  2.22  and  define  vl  to  be  the  local 
velocity  of  the  lattice  during  interdiffusion.  In  the  Kirkendall 
experiment,  after  some  annealing  time,  the  markers  plane  will  move 
in  the  positive  direction.  If  markers  are  placed  in  all  the  couple, 
those  in  the  diffusion  zone  will  move  in  a  different  extension.  In  the 
comers  there  will  not  be  movement  and  vl=  0,  whereas  in  the 

middle  (interface)  v/.  will  be  maximum.  From  his  analysis.  Darken 
derived  two  equations.  The  first  one: 

vl=[-Db  +  Da]  OXa/3x)  . (2.38) 

permits  to  obtain  the  Kirkendall  velocity  at  any  point.  For  example, 
at  the  comers  of  a  semi-infinite  diffusion  couple,  dKAldx  equals  zero, 
and  then  vi  =  0,  whereas  at  the  interface  the  local  gradient  will  be 

maximum  and  vi  will  be  maximum  too.  The  value  vi  can  be 

measured  experimentally  by  measuring  the  velocity  of  the  markers 

and  the  local  concentration  gradient  by  microprobe  analysis.  The 
other  Darken's  equation  can  be  expressed  as: 

D=XbDa  +  XaDb  .  (2.39) 

Again,  D  can  be  determined  experimentally  by  the  Boltzmann- 
Matano  method.  To  obtain  Da  and  Db  a  given  composition  one 
needs  to  determine  the  markers  velocity  and  the  interdiffusion  at 
the  position  in  the  couple  where  the  markers  reside.  From  the 
composition  profile  obtained,  for  instance,  by  microprobe  analysis 
the  variables  Xa>^b,  and  d^Aldx  can  be  determined  at  that  point. 
Then  by  solving  simultaneously  equation  2.38  and  2.39,  the  intrinsic 
diffusion  coefficients  can  be  calculated: 
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Da  =  D  +  Xa  VL  [3XA/axl 
Db  =  D  +  Xb  VL  [BXb/Bx] 


2.2.2 


(2.40) 


The  analysis  of  a  diffusion  process  in  multicomponent  systems 
(ternary  or  higher  order)  is  rather  complicated.  Even  in  the  most 
common  case  (temperature  and  pressure  constant,  no  external  force 
fields)  the  appearance  of  diffusion  flux  of  any  component  is  able  to 
cause  diffusion  redistribution  of  all  the  other  components.  As  to  the 
effort  needed,  in  general  the  amount  of  experimental  work  required 
increases  approximately  as  the  square  of  the  multiplicity  of  the 
system.  In  order  to  alleviate  the  complexity  of  the  analysis,  "the 
experimental  approach  of  ternary  systems  must  be  accompanied  by 
an  exhaustive  theoretical  investigation  of  all  aspects  of  the  problem" 
[64].  In  this  section  a  brief  review  of  the  formal  elements  of  the 
ternary  diffusion  theory  in  single  phase  substitutional  systems  is 
presented.  For  the  sake  of  simplicity  during  the  analysis  possible 
effects  such  as  volume  changes  of  mixing,  lateral  dimensional 
changes  and  pore  formation  will  be  neglected. 


2.2.2. 1  The  Diffusivitv  Matrix 


At  present  the  most  common  method  to  define  the  total  set  of 
mutual  diffusion  coefficients  in  multicomponent  systems  is  that  of 
Matano-Kirkaldy  [64).  A  brief  explanation  of  its  fundaments  follows. 


As  in  the  binary  case,  the  flux  equations  may  be  written  in  terms  of 
either  the  chemical  potential  gradient  (cal/mole)  or  the  concentration 
gradient  (mole/cm3)  in  the  Kirkendall  (lattice  fixed)  frame  of 
reference.  Generally,  the  flow  equation  (Pick's  First  Law)  for  this  case 
is: 

c 

Ji=SlDij]Vcj  .  (2.41) 

j-2 

where  [Dij]  is  known  as  the  matrix  of  diffusivities.  Note  that  yruns 
from  2  to  c  in  this  equation.  In  a  ternary  system  there  are  two 
independent  composition  variables.  Therefore,  a  diffusion  couple 
experiment  requires  the  determination  of  two  independent 
composition  profiles  to  characterize  a  couple.  Since  to  a  good 


approximation: 

Cl  +  C2  +  C3  =  c  (moles/cm^)=  constant  .  (2.42) 

and  its  derivative: 

Vci  +  Vc2  +  Vc3  =  0  .  (2.43) 

then,  the  dependent  variable: 

Vc3  =  -(Vci  +  Vc2  )  .  (2.44) 

Regarding  with  the  index  i  in  the  equation  2.41,  it  runs  only  over  two 
components  because  the  3  fluxes  are  not  independents: 

Jl  +  J2  +  J3  =  0  .  (2.45) 


From  Gibb’s  phase  rule,  ternary  systems  offer  an  extra  degree 
of  freedom  as  compared  to  binary  couples.  This  fact  means  that  the 
sequence  of  compositions  (so-called  composition  path)  in  ternary 
systems  may  be  different  from  couple  to  couple  in  the  same  system. 
Such  composition  path  is  time  invariant  for  semi-infinite  couples  and 
is  represented  by  a  continuous  curve  on  the  Gibbs  triangle  ternary 
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phase  diagram  representation  that  connects  the  endpoint 
compositions  of  the  alloys  making  up  the  diffusion  couple  (Figure 
2.23). 


When  the  dependent  variable  is  component  3,  equation  2.41 

can  be  expressed  as: 

Jl  =  -  Dii^3ci/9x  ■  Di2^^C2/3x 


.  (2.46) 

J2  =  -  D2i^3ci/9x  *  D22^^C2/9x 

One  problem  associated  with  the  definition  of  the  diffusion 

coefficient  in  a  ternary  system,  however,  is  that  it  depends  upon  the 

choices  made  for  the  dependent  flux  and  concentration  gradient;  so, 

the  same  equation  2.41  can  be  rewritten  as  (using  short 

nomenclature): 

Jl  =  -  Dn^Vci  -  Dij^Vcs 


J2=-D2,2Vci-D23'Vc3 


(2.47) 


when  component  2  is  the  dependent  variable.  Substituting  2.44  in 
2.47: 

Jl  =  -  Dn^Vci  -  D,3^(-Vci  -  Vc2  ) 

J2  =  -  ^21^  Vci  -  D23^(-Vci  -  VC2  ) 

then. 


Jl  =  -  [D„2-  D,3^]  Vci  -  DnH'Vci  -  Vc2  ) 

J2  =  -  [D2,^-  D23']  Vci  -  D23^(-Vci  -  VC2  ) 

comparing  with  2.46: 

D„'=  [D„2-D,3^] 


D,2^=-D 


D2,"=  [D2,"-  D23"] 
D22^=  -  D23^ 


(2.48) 


(2.49) 


These  relationships  among  the  coefficients  of  diffusion  lack  physical 
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meaning,  and  a  result  such  as  D;;>D22  may  not  necessarily  imply  that 
component  1  moves  faster  than  component  2  in  the  system. 

In  equation  2.46,  the  on-diagonal  element  (<=y>,  known  also  as 
the  mutual  diffusion  coefficients,  reflect  how  the  gradient  of 
component  i  influences  the  flux  of  the  same  component.  The  so-called 
off  diagonal  terms  (or  indirect  coefficients)  report  how  the  gradient 
of  the  other  components  affect  the  flux  of  the  component  / . 
Frequently,  the  off-diagonal  terms  are  much  smaller  than  the  on- 
diagonal  terms. 

2.2.2.2  Experimental  Determination  of  fDiil.  The  infinite  Diffu.sion 
couple. 

As  in  the  binary  case,  the  infinite  diffusion  couple  is  one  of  the 

most  useful  techniques  to  obtain  solid-state  ternay  diffusion  data 

because  the  associate  boundary  conditions  leads  to  a  particular 

rather  simply  integrals  of  the  diffusion  equation.  For  instance, 

applying  local  conservation  of  mass  on  equation  2.46  it  is  possible  to 

derive  the  equivalent  of  Fick's  Second  Law: 

5ci/at=  3(Dii^3ci/ax)/ax  +  ^(Dl2’^C2/^x)/^x 

.  (2.50) 

3c2/at=  5(D2i^3ci/ax)/ax  +  ^(I>22^^c2/ax)/ax 

subject  to  the  boundary  and  initial  conditions: 

ci(±  X,  0)  =  Ci(±  oo,  t)  =  Ci±  .  (2.5 1 ) 


As  in  the  binary  system,  in  this  case  it  is  possible  to  make  the  lamda- 
substitution  (equation  2.17),  and  convert  the  partial  differential 


equation  2.50  into  a  total  differential  equation,  which  after  the 

Matano-Boltzmann  solution  yields: 

c,a) 

-  (l/2t)  J  X  dci  =  D„'dci/dX  I  +  D,2^dc2/dX  I  .  (2.52a) 

Cl-  Ci(X)  C2(X) 

similarly, 

C2W 

-  (l/2t)  J  dci  =  D2iMci/dX  I  +  D22Mc2/dX  |  .  (2.52b) 

C2-  Ci(X)  C2(X) 

The  left  side  of  equations  2.52  can  be  empirically  evaluated  with  the 
composition  profile  of  each  component  (Matano  Integral).  This  is 
graphically  illustrated  in  Figure  2.24.  Then  for  a  particular 
composition,  the  information  needed  to  evaluate  equation  2.52  is 
fXdc^  obtained  by  a  numerical  technique,  and  the  slope  of  the  curve 
at  that  composition  point,  which  can  be  readily  obtained.  Thus,  at  any 
composition  there  are  two  simultaneous  linear  equations  with  a  set 
of  4  unknowns:  D/y,  D/2,  D21.  D22-  To  make  the  system  determinate, 
therefore,  it  is  necessary  to  produce  another  independent  diffusion 
couple  and  perform  the  same  analysis  on  its  composition  profile. 

Consider  Figure  2.23.  If  the  diffusion  couples  have  been  chosen 
so  that  their  compositions  paths  on  the  ternary  isotherm  cross  at 
coordinates  (cy,C2),  we  see  that  there  are  4  equations  and  4  unknowns 
soluble  only  at  that  concentration  {C1C2)  which  is  common  to  both 
couples;  the  values  of  the  set  of  form  Dy/j,  being  functions  only  of 
compositions  but  not  the  gradients,  will  be  the  same  for  that 
composition  that  appears  in  both  couples.  Thus,  the  bottom  line  is 
that  the  matrix  /D,y7  can  be  determined  by  the  Matano-like  analysis 
in  two  diffusion  couples  at  their  crossing  points,  over  which  no 
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control  exists.  One  problem  associated  is  that  the  crossing  point 
composition  is  usually  one  end  of  the  composition  path  which  makes 
no  easy  the  evaluation  of  the  integral  and  slope,  diminishing  the 
precision  of  the  result.  Then,  the  general  strategy  employed  in 
experimental  studies  generates  a  whole  collection  of  diffusion 
couples  that  survey  the  ternary  region  to  be  studied.  After 
determining  the  resultant  composition  profiles  and  their 
corresponding  composition  paths,  the  integrals  and  derivatives  are 
evaluated  at  all  crossing  points  that  happen  to  develop  between  pairs 
of  couples.  It  should  be  pointed  out  that  a  given  couple  could  have  a 
lot  of  intersections  with  other  couples.  The  results  of  this  kind  of 
study  are  usually  reported  as  a  counter  map  of  each  of  the  diffusion 
coefficients:  Du,  Du,  D21,  D22  as  a  function  of  composition  on  the 
Gibbs  triangle  (Figure  2.25)  [65].  These  values  are  far  from  the  real 
physical  meaning,  so  that  if  D// >£>72  that  does  not  mean  necessarily 
that  component  1  moves  faster  than  component  2.  The  interdiffusion 
studies  in  Cu-Ag-Au  system  by  Ziebold  and  Oglivie  are  a  typical 
application  of  these  strategy  [65]. 

The  amount  of  experimental  work  required  to  survey  the  2- 
dimensional  surface  of  the  Gibbs’  triangle  for  each  Dij  can  be  reduced 
drastically  by  the  use  of  special  procedures  and  theoretical 
knowledge  about  the  diffusivities.  Special  procedures  for  evaluating 
the  3  independent  mutual  coefficients  have  been  described  by  [66- 
68]. 

Some  important  experimental  tendencies  in  typical  ternary 
studies  indicate  that  the  maximum  in  D  occurs  usually  close  to  the 
corresponding  value  of  the  minimum  melting  point.  Another 
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GOLD 


Figure  2.25  a)  The  direct  coefficient  of  diffusion  for  Cu  Dn^ 
in  units  10’*  2  cm^/sec  at  725  °C.  Filled  circles 
indicate  intersections  of  diffusion  paths,  and 
open  squares  indicate  extreme  points  in  the 
composition  profile  [65]. 
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Figure  2.25 


(Continue)  b)  The  direct  coefficient  of  diffusion 
for  Ag  D22^  in  units  10‘12  cm^/scc  at  725  ®C. 
Filled  circles  indicate  intersections  of  diffusion 
paths,  and  open  squares  indicate  extreme  points 
in  the  composition  profile  [65]. 
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interesting  observation  is  that  to  calculate  D/yfrom  extreme  values 
(maximum  or  minimum)  in  the  composition  profiles,  we  do  not  need 
a  crossing  point  since  the  Vaj  equals  to  zero  at  that  particular 
concentration  reducing  the  equations  2.46  to  a  system  of  2  equations 
and  2  unknowns: 

Jl  =  -  Du^3ci/ax  -  D,2^5c2/ax  =  -  DM^^Cl/^x 

.  (2.53) 

j2  =  -  D2i^aci/ax  -  D22^ac2/ax  =  -  D2i^aci/ax 


2. 2.2.4  Diffusion  Penetration  Tendencies  and  Composition  Path 

In  1974,  DeHoff  et.  al.  [69]  published  a  new  viewpoint  on  the 
analysis  of  the  diffusion  behavior  in  ternary  systems  where  [Dij]  is 
replaced  by  the  relative  penetration  tendencies  of  the  3  species. 
Through  this  approach,  the  calculation  of  a  complete  pattern  of 
composition  paths  and  Kirkendall  shifts  on  the  Gibbs  triangle  for  the 
relative  penetration  tendencies  of  the  components  is  relatively 
simple. 

During  their  analysis,  they  made  3  assumptions.  First,  the 
interpenetrations  of  the  species  occur  on  a  rigid  lattice  so  that 
appreciable  excess  and  deficit  of  atoms  can  occur  on  one  side  and  the 
other.  To  do  this,  the  normal  continuous  feedback  of  vacancies  was 
changed  for  a  one-step  feedback  that  permitted  the  readjustment  of 
the  lattice.  This  was  recognized  as  physically  unreal,  but  essential  for 
the  calculation  of  the  Kirkendall  shift  for  the  system.  The  third 
element  of  the  approach  was  approximating  the  concentration- 
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penetration  curves  for  each  component  as  straight  line  with  the  rigid 
lattice  as  the  coordinate  system.  In  this  case,  when  the  depth  of 
penetration  of  an  element  depended  on  composition,  the 
concentration-penetration  curve  was  approximated  by  two  straight 
line  segments. 

Consider  a  diffusion  couple  made  up  by  alloys  A-50at%B  (P 
side)  and  B-50at%C  (Q  side)  in  a  single  ternary  system  A-B-C  whose 
relative  penetration  depths  are  independent  of  composition  and 
already  known:  A>B>C.  Figure  2.26a  shows  the  linear  composition 
profiles  (in  atoms/cm^)  for  this  couple  in  a  rigid  lattice  before  and 
after  the  diffusion  treatment.  In  accordance  with  the  relative  depths 
of  penetration,  component  A  draws  deeply  from  the  P  side  of  the 
couple  and  penetrates  deeply  in  the  Q  side  of  the  couple.  Regardless 
of  its  mobility,  component  B  does  not  move  since  no  concentration 
gradient  for  it  is  present.  Component  C  slightly  penetrates  into  the  P 
side  and  draws  from  a  shallow  depth  on  the  Q  side.  The  total 
concentration  of  atoms  in  the  unshifted  lattice  (Figure  2.26b), 
obtained  by  adding  the  number  of  atoms  of  the  3  species  point  by 
point,  clearly  shews  a  deficiency  of  atoms/cm^  in  the  P  side  and  an 
excess  of  them  in  the  Q  side.  Since  co  is  considered  to  be  independent 
of  composition,  it  is  necessary  to  generate  a  number  of  vacancies 
proportional  to  the  area  A  j  in  the  right  side,  which  in  turn  will 
increase  the  thickness  of  the  Q  side  and  decrease  the  P  side  of  the 
couple.  This  unbalance  will  produce  a  Kirkendall  shift  (translation  of 
the  original  weld  interface)  toward  the  left  side  of  the  couple,  shown 
as  Xk  in  Figure  2.26c.  Since  all  the  positions  in  the  diffusion  zone  are 
shifted  by  this  process,  the  original  coordinate  system  x<>  (rigid 


68 


Figure  2.26  a)  Model  for  a  ternary  diffusion  couple  with  independent 

concentration  penetration  tendencies;  b)  total  concentration  curve; 
c)  computation  of  Kirkendall  shift;  d)  prediction  of  Kirkendall 
shift:  e)  atom-fraction  penetration  profiles  for  the  unrelaxed 
lattice:  and  f)  computed  composition  path  (69]. 
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lattice)  is  transformed  to  a  new  coordinate  system  x  which  is  all 
moved  to  the  left.  The  magnitude  of  this  change,  however,  is 
different  in  different  positions.  For  example,  at  xy^this  magnitude  is 
proportional  to  the  volume  change  accumulated  to  the  left  of  this 
point  by  the  annihilation  of  vacancies  (area  A3),  which  must  be  equal 
to  the  to  the  volume  occupied  by  these  vacant  sites  prior  to  their 
annihilation  (area  A 4).  By  equating  these  two  areas  at  each  point  of 
the  affected  lattice,  it  is  possible  estimate  the  magnitude  of  the 
lattice  displacement  of  the  original  system.  Ax.  Note  that  the 
Kirkendall  shift  is  maximum  at  the  weld  interface  (Figure  2.26d). 

The  atomic  fractions  of  Xa,Xb.  and  Xc  are  shown  in  Figure 
2.26e.  They  are  derived  at  any  point  of  the  unshifted  system  by 

dividing  the  corresponding  c/  by  the  respective  Ct.  This  curves  would 
be  the  experimental  result  given  by  the  microprobe  analysis. 
Observe  how  the  profile  corresponding  to  the  component  B  varies  up 
and  down  as  it  crosses  the  diffusion  zone  suggesting  a  movement  of 
atoms  B  from  the  right  side  of  the  couple  in  a  possible  uphill 
diffusion  way.  However,  as  mentioned  earlier  there  is  not  effective 
flow  of  atoms  B  at  all  because  of  the  absence  of  a  B  concentration 
gradient,  and  the  variation  of  its  composition  is  due  only  to  the 

imbalance  in  the  distribution  of  A  and  C  atoms.  Finally  Figure  2.26f 

shows  a  plot  of  the  composition  path  and  Kirkendall  shift  that  would 
be  observed  for  this  kind  of  system.  There  are  some  interesting 

observations  that  should  be  noted  in  this  plot.  For  example,  the  two 
ends  of  the  path  point  directly  toward  the  corner  of  the  fastest 
component  A,  and  the  central  part  of  the  path  tends  to  point  toward 
the  corner  of  the  slowest  component  C.  This  causes  a  considerable 
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deviation  of  the  composition  path  from  the  straight  line  joining  the 
endpoint  compositions  (concentration  tie-line)  with  the  path  curve 
crossing  at  half  way  of  this  tie-line.  The  intersecting  point  is  called 
the  cross  over-point  and  is  important  in  the  determination  of  the 
Kirkendall  shift. 

When  the  mobility  of  the  atomic  species  in  the  couple  varies 
with  composition,  the  depths  of  penetration  and  draw  for  each 
component  are  different  on  the  2  sides  of  the  couple.  This  case  is 
illustrated  in  Figure  2.27  [69].  The  only  difference  with  the  previous 
case  is  that  the  penetration  and  draw  depths  are  larger  in  the  left 
side  of  the  couple  (richest  in  the  fastest  component  A )  and  more 
shallow  in  the  other  side.  The  determination  of  the  concentration  at 
the  welded  interface  is  given  by  the  fact  that  the  triangular  areas  on 
both  sides  of  the  couple  must  be  equal  for  each  component 
(conservation  of  atoms).  A  similar  analysis  to  that  given  above  would 
explain  the  plots  for  this  case.  Nevertheless,  some  interesting 
deviations  in  the  composition  path  for  this  case  must  be  mentioned. 
One  is  the  displacement  of  the  crossover  point  toward  the  P  side  of 
the  couple.  An  important  characteristic  of  the  crossover  point  is  that 
its  composition  correspond  to  that  of  the  original  weld  interface.  The 
other  difference  is  that  the  penetration  curves  are  skewed,  which 
causes  an  enlargement  and  a  shortening  in  the  the  length  of  the 
segments  at  the  end  of  the  composition  path  in  the  P  and  Q  side, 
respectively. 

This  approach  was  successfully  applied  by  the  authors  to 
generate  simple  penetration  tendency  models  that  reproduce 


Figure  2.27  a)  Model  for  a  ternary  diffusion  couple  with  concentration  dependent 
penetration  tendencies;  b)  toul  concentration  curve;  c)  computation  of 
Kirkendall  shift;  d)  prediction  of  Kirkendall  shift;  e)  atom-fraction 
penetration  profiles  for  the  unrelaxed  lattice;  and  0  computed 
composition  path.  Note  the  ondulation  in  component  B  (e)  and  the 
displacement  of  the  crossover  point  on  the  path  (0.  [69] 
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experimental  composition  paths  and  Kirkendall  shift  patterns  in 
different  ternary  systems,  such  as  Cu-Ag-Au,  Cu-Ni-Zn  and  Cu-Ni-Zn. 

2.2.2.3  Zero  Flux  Plane 

In  1979,  Dayananda  and  Kim  [70]  described  a  different  way  to 
calculate  interdiffusion  fluxes  in  binary  and  multicomponent  systems 
directly  from  composition  profiles  of  isothermal  diffusion  couples  by 
an  analysis  directly  derivable  from  the  continuity  equation  without 
the  need  for  a  prior  knowledge  of  interdiffusion  coefficients,  the  only 
assumption  being  constant  molar  volume.  A  major  outcome  of  their 
interdiffusivity  studies  on  the  Cu'Ni-Zn  system  was  the  identification 
of  ’zero  flux  plane'  for  each  component  within  the  diffusion  zone  of 
the  diffusion  couples.  They  defined  a  ’zero  flux  plane’  as  the  plane 
where  the  interdiffusion  flux  of  a  component  goes  to  zero  and  on 
either  side  of  the  plane  occurs  a  change  or  reversal  in  the  direction  of 
the  interdiffusion  flux  of  the  component. 


CHAPTERS 

MATERIALS  AND  EXPERIMENTAL  PROCEDURES 

In  this  chapter  a  description  of  materials  and  experimental 
techniques  employed  throughout  this  research  is  given.  A  flow  chart 
that  represents  the  experimental  procedure  for  preparation  and 
analysis  of  diffusion  couples  is  given  in  Figure  3.1.  The  experimental 
procedure  consists  of  alloy  preparation  by  arc-melting, 
homogenization,  diffusion  couple  preparation,  interdiffusion  heat 
treatment,  analytical  analysis  using  microprobe,  and  computerized 
analysis  of  the  composition  profiles.  These  steps  are  explained  in  the 
following  sections. 

3.1  Materials 

The  alloys  used  in  this  study  were  prepared  by  arc-melting  of 
high  purity  materials  (99,99%  Nb  and  Ti,  99.999%  Al).  The  pure 
metals  were  in  the  form  of  a  5-mm-thick  niobium  sheet  obtained 
from  Teledyne-Wah-Chang,  a  6.35-mm-diameter  titanium  rod  from 
Aesar,  and  a  38-mm-diameter  aluminum  rod  from  United  Mineral 
Corporation.  Table  3.1  gives  the  ingot  analysis  of  the  pure  metals. 
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Figure  3.1-  A  flow  chan  showing  ihc  cxpcrimenial  procedure  for 
the  preparation  and  analysis  of  the  diffusion  couples. 


3.2  Alloy  Preparation 


The  binary  and  ternary  alloys  were  prepared  by  arc  melting 
weighed  amounts  of  pure  metals  on  a  water-cooled-jacket  copper 
stage  under  an  oxygen-gettered  argon  atmosphere.  Before  melting,  a 
vacuum  of  100  mtorr  in  the  chamber  was  attained  and  then  argon, 
previously  dried  in  a  gettering  furnace  was  allowed  into  the  system 
up  to  a  pressure  of  -15  KPa.  The  purging  operation  was  repeated  at 
least  4  times  before  proceeding  to  melt.  The  specimens  were  melted 
under  an  Ar  pressure  of  -35  KPa.  Each  specimen  was  melted  at  least 
4  times  to  assure  compositional  homogeneity  of  the  buttons.  The 
buttons  were  weighed  after  each  melting  operation  so  that  A1  losses 
could  be  monitored.  Weight  losses  were  detected  only  in  buttons  of 
the  binary  Nb-Al  alloy  and  accounted  for  approximately  25%  of  the 
initial  A1  weight.  In  order  to  obtain  the  desired  composition  in  the 
alloy,  an  excess  of  pure  Al,  equivalent  to  these  losses,  was  added  to 
the  samples  before  melting. 

After  metallography  preparation,  the  alloys  were  analyzed  by 
microprobe  analysis,  whose  results  proved  that  the  Al-loss 
estimation  was  correct.  The  resultant  compositions  from  the 
microprobe  analysis  are  shown  in  Table  3.2.  In  order  to  check  their 
accuracy,  some  of  the  alloys  were  sent  out  for  wet  chemical  analysis. 
Table  3.3  compares  the  results  obtained  by  both  techniques  for  the 
binary  alloy.  A  good  agreement  of  results  is  observed. 


Table  3.2.  Alloy  composilions 
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3.3  Preparation  of  Diffusion  Couples 

The  arc  melted  buttons  in  the  form  of  spherical  cap-like  shape 
of  approximately  3  gram  weight  and  8  mm  height,  were  studied 
using  optical  and  electron  microscopy.  Microprobe  analysis  showed 
macrosegregation  and  significant  microsegregation  in  the 
microstructure.  Therefore,  a  homogenization  of  the  alloys  was 
required.  The  buttons  were  plastically  deformed  to  a  disk-like  shape, 
5 -mm  thick,  at  room  temperature  using  a  press.  The  new  shape 
permitted  to  obtain  flat  surfaces  and  the  plastic  deformation 
expedited  the  subsequent  homogenization  through  breaking  down  of 
the  dendrite  arms.  All  the  heat  treatments  performed  in  this 
research,  i.e.  the  homogenization  and  the  diffusion  annealings,  were 
done  in  a  Minivac  Furnace  Series  2400  Vacuum  Bell  Jar  ,  which  is  an 
ultra-high  temperature  vertical  cold  wall  vacuum  furnace  capable  of 
temperatures  up  to  2000  ®C.  A  schematic  representation  of  the 
furnace  assembly  employed  for  all  the  heat  treatments  is  illustrated 
in  Figure  3.2  The  pumping  system,  consisting  of  mechanical  and 
diffusion  pumps  and  furnished  with  ionization  gauges,  is  able  to 
produce  a  10'^  torr  vacuum  in  the  bell  jar.  To  prevent  any 
penetration  of  oxygen  into  the  alloys  during  heat  treatments,  the 
samples  were  wrapped  in  tantalum  foil  (100  pm  thick).  The  binary 
Nb-Al  alloy  (po)  proved  to  be  also  an  excellent  oxygen  getter.  Before 
heating  up  the  samples,  a  vacuum  of  lE-5  was  obtained  in  the  bell 
Jar.  The  conditions  for  the  homogenization  treatment  depended  on 
the  chemical  composition  of  the  samples.  In  the  case  of  binary  alloy 
Po,  the  treatment  was  done  at  1400  ®C  for  6  hours.  For  all  other 


Figure  3.2  Furnace  assembly  used  for  all  beat  (reatments 
in  this  research. 


alloys,  the  homogenization  was  conducted  at  the  same  temperature 
but  for  only  4  hours.  The  lower  time  is  due  to  the  presence  of 
titanium  which  reduces  the  melting  point,  and  subsequently 

increases  the  rate  of  diffusivity  and  the  homogenization  process.  The 
system  was  heated  up  to  950  °C  for  30  minutes,  and  then  brought  to 
the  homogenization  temperature  of  1400  °C  in  5  minutes.  After 
annealing,  the  samples  were  cooled  down  in  the  furnace.  Optical 
microscopy  and  backscattered  electron  analysis  were  performed  on 
the  heat  treated  samples  to  ensure  that  the  microsegregation  was 

removed  and  specimens  were  homogenized.  Since  an 
environmentally  affected  zone,  which  is  depleted  in  A1  and  has  high 
oxygen  content,  develops  in  the  vicinity  of  the  surface  of  the 
homogenized  specimens,  the  surface  layer  was  removed  and 

rectangular  specimens  of  (1  x  I  x  0.5  cm^)  were  cut  from  the 

homogenized  buttons.  Samples  were  prepared  metallographically  for 
the  subsequent  diffusion  bonding  using  a  fixture  that  permitted  to 
obtain  parallel  surface. 

As  mentioned  in  section  2.3,  one  of  the  simplest  techniques  to 
assemble  solid  state  diffusion  couples  is  the  window-frame  technique 
[71].  In  this  method,  the  two  members  pf  the  couple  are  polished  flat 
and  mechanically  clamped  together.  In  some  cases,  like  in  this  study, 
a  short  heat  treatment  at  a  temperature  below  the  interdiffusion 
heat  temperature  is  used  to  facilitate  the  diffusion  bonding.  The  use 
of  suitable  markers  to  delineate  the  original  bonded  interface  is 
important.  Their  size  should  be  small  enough  to  avoid  interfering 
with  the  diffusion  bonding  process  but  big  enough  to  be  detected  by 
the  SEM.  They  also  should  be  chemically  inert  with  the  components 
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of  the  couple  to  avoid  their  decomposition  during  the  heat 
treatments.  For  this  research,  yttria  (Y2O3)  particles,  1  pm  in  size, 
were  employed. 

The  diffusion  couples  for  this  investigation  were  prepared 
using  a  Centorr  Series  600  Hot-Press  Furnace.  The  steel  chamber  has 
a  workable  space  of  76  mm  diameter  x  127  mm  high  and  is  water 
cooled  by  a  flow  of  5  GPM  at  50  psi  and  20  °C.  This  system  uses  an 
induction  coil  as  a  heater  element  to  achieve  temperatures  up  to 
1500  “C  in  vacuum  and  neutral  atmospheres.  It  is  equipped  with  a 
mechanical  pump  that  permits  to  achieve  a  vacuum  in  the  order  of 
lE-2  torr.  For  inert  gas  operations,  inlet  valves,  relief  valves  and  a 
combination  pressure-vacuum  gauge  are  provided.  To  operate  with 
argon  or  any  other  inert  gas  atmosphere,  the  chamber  can  be 
evacuated  by  means  of  the  pumping  system  and  then  backfilled  with 
the  desired  inert  gas.  The  specimens  to  form  the  diffusion  couple 
were  compressed  together  by  means  of  a  hydraulic  press  with  a  steel 
ram,  25.4  mm  in  diameter,  coming  down  from  the  top  of  the 
chamber.  The  temperature  measurements  were  carried  out  by  a  fast- 
respond  optical  pyrometer  with  adjusted  chromatic  compensation 
and  a  precision  of  ±  10  °C. 

Two  kinds  of  diffusion  couples  were  built  for  this  research.  For 

the  binary  diffusion  study  a  two-part  couple  was  used,  whereas  for 

the  ternary  diffusion  case  a  sandwich-like  couple  was  made.  Small 

drops  of  a  dilute  solution  of  the  fine  chopped  yttria  particles  in 
acetone  were  deposited  as  markers  over  the  polished  surface  of  the 
specimens.  The  couple  was  compressed  together  inside  a  graphite 

die.  The  3-piece  die  consists  of  a  stage  where  the  sample  to  be 
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compressed  rests,  a  ram  through  which  the  load  is  applied,  and  a 
case  to  hold  the  whole  system.  Before  setting  the  sample  in  the 
graphite  die,  a  very  thin  film  of  boron  nitride,  a  very  stable  material 
at  the  employed  bonding  temperatures,  was  sprayed  over  the  faces 
in  contact  with  the  specimen  so  that  the  diffusivity  of  the  carbon  into 
the  sample  was  prevented.  In  addition,  a  lOO-pm-thick  tantalum  foil 
was  placed  between  the  boron  nitride  layer  and  the  sample.  The 
whole  system  was  then  set  into  the  induction  coil  in  the  chamber,  as 
shown  in  Figure  3.3.  As  in  the  arc  melting  operation,  the  chamber 
was  purged  5  times  with  gettered  Ar  before  heating  up  the  system  to 
the  bonding  temperature. 

The  diffusion  couples  investigated  in  this  research  are  listed  in 
Table  3.3.  The  bonding  temperature  was  1030  ±  10  ®C  for  the  binary 
diffusion  system  Nb(p)/|30  and  930  ±  10  ®C  for  the  couples  including 
ternary  alloys,  except  the  couple  made  with  pure  Ti.  Once  at  the 
bonding  temperature,  a  stress  of  approximately  2  MPa  was  applied 
to  the  system  and  held  for  the  length  of  the  diffusion  bonding.  After 
cooling  down,  the  bonded  diffusion  couple  was  cut  into  two  parts.  A 
section  perpendicular  to  the  bonded  interface  was  metallographically 
prepared  for  examining  the  quality  of  the  bonded  interface.  The 
characterization  of  the  bonded  interface  prior  to  the  interdiffusion 
was  done  using  optical  microscopy,  scanning  electron  microscopy  and 
microprobe  analysis.  A  solution  of  Kroll’s  reagent  (10  ml  HF,  30  ml 
H2NO3,  50  ml  H2O)  was  used  as  an  etchant  for  revealing  the 
microstructure  of  the  ternary  Nb-Al-Ti  alloys. 


84 


LOAD 

4' 


GRAPHITE 

DIE 


■Boron  Nitride - 
Tontalum  Foil 
Layers  Set 


Figure  3.3 


Experimental  set-up  for  hot-pressing  bonding  of 
diffusion  couples. 
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3.4  Interdiffusion  Heat  Treatments 

The  diffusion  couples  were  wrapped  in  tantalum  foil,  put  upon 
the  alumina  stage  in  the  vacuum  furnace,  as  shown  in  Figure  3.2,  and 
were  surrounded  by  the  Nb-Al  binary  alloy  as  a  getter.  A  vacuum  of 
lE-5  torr  was  set  in  the  bell  jar  before  heating  up  to  the 
interdiffusion  temperature.  Again,  similar  to  the  homogenization 
process,  a  short  heating  time  («5  minutes)  from  950  ®C  to  the 

treatment  temperature  was  registered.  A  close  control  of 
temperature  (±5  °C)  was  achieved  by  using  a  type  C  (W-5%Re) 
thermocouple.  Since  the  vacuum  bell  jar  furnace  lacks  a  quenching 
capability,  the  samples  were  cooled  inside  the  furnace.  The  cooling 
rate,  however,  was  fast  enough  to  avoid  further  diffusion  upon 
cooling.  For  example,  it  took  less  than  90  seconds  to  go  from  1550®  to 
950  ®C,  which  is  a  nr  'ligible  time  compared  to  the  duration  of  the 
interdiffusion  treatments  employed  in  this  research. 

To  evaluate  the  diffusivity  of  aluminum  in  the  binary  Nb-Al 
alloys,  the  couples  (#1  to  #5)  were  annealed  for  various  lenght  of 
time  in  the  temperature  range  of  1350®- 1550  ®C.  The  duration  of  the 
heat  treatment  was  limited  by  a  lower  bound  sufficient  to  produce  a 
diffusivity  zone  of  at  least  25  pm  at  the  lowest  treatment 
temperature,  and  an  upper  bound  to  avoid  alumina  formation  in  the 
binary  Nb-Al  alloy  at  higher  temperatures.  Table  3.4  provides  the 
details  regarding  the  binary  interdiffusion  treatments  conducted. 

With  respect  to  the  ternary  system,  diffusion  treatments  of  the 
sandwich-like  diffusion  couples  were  conducted  at  two  different 
temperatures  (1200®  and  1400  ®C).  First,  diffusion  couples  #6  and  #7 
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were  treated  at  1400  ®C.  Since  the  resultant  penetration  zone  was 
very  large,  the  treatment  temperature  was  decreased  to  1200  °C.  The 
interdiffusion  time  for  the  ternary  alloys  was  held  constant  (6  hours) 
in  order  to  compare  the  penetration  tendencies  in  these  couples. 

3.5  Analytical  Microscopy 

The  diffusion  couples  used  in  this  study  as  well  as  the 
materials  employed  during  their  construction  were  chemically 
analyzed  and  structurally  characterized  by  Microprobe  Analysis 
(EPMA),  and  Scanning  Electron  Microscopy  (SEM)  (Backscattering 
Electron  BSE  and  Secondary  Electron  SE  modes),  respectively.  These 
analyses  were  performed  in  the  diffusion  couples  before  and  after 
the  diffusion  treatment  in  order  to  check  the  extent  of  diffusion  and 
possible  phase  transformation  during  the  bonding  process,  and  to 
obtain  the  respective  concentration  profiles  for  the  diffusivity 
analysis.  Concentration  profiles  were  obtained  using  a  JEOL 
Superprobe  733  electron  micrcprobe.  This  system  is  equipped  with  a 
tungsten  filament  electron  gun  and  four  wavelength  dispersive 
spectrometers  (WDS).  A  15-kV,  l-pm  electron  beam  with  a  current  of 
100  nanoamp  was  employed  in  this  research.  The  quantitative 
interpretation  of  the  x-ray  signal  is  obtained  after  applying  the  ZAP 
technique,  which  is  based  on  standard  element  intensities  and 
corrects  the  measured  intensities  for  experimental  effects  such  as 
atomic  number  Z,  absorption  A  and  fluorescence  F  [72].  The 
compositions  of  Nb,  Al,  and  Ti  were  calculated  by  monitoring  and 
measuring  the  intensities  of  AlKa  (E5=1.487  keV,  \=  8.339  A),  TiKa  (E= 
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4.511  keV,  X=  2.749  A),  and  NbLa  (£=2.166  keV,  X=  5.724  A)  x-ray 
signals.  The  Nb-  and  Ti-  x-ray  radiations  were  detected  using  a 
system  of  PET  crystal/Xe  gas  filled  detector.  The  PET  crystal  has  a 
Bragg’s  equation's  characteristic  2d  value  of  8.742  A  and  is  usually 
used  for  relatively  long  wavelength  radiations.  On  the  other  hand, 
the  A1  x-ray  signals  were  detected  using  a  TAP  (Thalium  Acid 
Phthalate)  analyzing  crystal  with  a  characteristic  2d  value  of  25.757 
A,  and  a  PTlO-gas-flow  (90%  Ar-10%  CH4)  proportional  counter.  This 
system  is  commonly  used  for  detecting  the  x-ray  energies  of  light 
elements.  The  x-ray  counting  process  during  the  microprobe  analysis 
was  statistically  controlled.  A  composition  standard  deviation  of  0.1 
at%  Al,  and  a  200-second  maximum  counting  time  were  set  to  obtain 
all  the  composition  profiles  in  this  investigation.  The  concentration 
profiles  were  determined  by  a  point-to-point  counting  technique.  For 
the  binary  diffusion  case  the  measurements  were  done  using  1  pm 
steps.  Larger  steps  (5-10  pm)  were  required  for  those  ternary 
diffusion  couples  (1400  °C  treatments)  which  showed  large 
penetration  zones. 

Microstructural  analysis  of  the  samples  were  performed  in  a 
SEM  Jeol  model  103.  BSE  analysis  was  conducted  in  samples  where 
chemical  composition  issues  were  wanted,  whereas  SE  mode  was 
used  when  the  morphology  of  the  structure  was  being  analyzed.  This 
microscope  has  the  capability  to  detect  elements  with  Z>5  for 
qualitative  analysis. 
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3.6  Calculation  Method 

The  composition  profile  data  obtained  from  the  binary 
diffusion  couples  were  then  introduced  to  a  Fortran-77  source  code 
[8]  that  includes  several  subroutines.  The  first  of  them  was  used  to 
fit  a  smooth  composition  curve  to  the  experimental  composition 
profiles.  Then,  a  different  subprogram  determined  the  position 
where  the  Matano  Interface  was  located  in  the  fitted  curve,  and 
finally,  the  last  subroutine  calculated  the  concentration  gradient,  the 
interdiffusion  flux,  and  the  coefficients  of  interdiffusion  on  the  basis 
of  the  fitted  concentration  curve  and  the  Matano  plane  location.  Each 
of  the  subprograms  makes  use  of  routines  from  the  International 
Mathematics  and  Statistics  Library  (IMSL).  For  the  fitting  of 
experimental  data,  a  modified  cubic  spline  routine  (ICSVKU)  was 
employed.  The  modifications  to  such  routine  include  the  introduction 
in  each  side  of  the  concentration  profile  of  a  linear  zone,  which  spline 
cubic  equations  can  not  account  for,  and  a  parabolic  zone,  which 
smooths  the  connection  between  the  linear  and  cubic  curves.  Details 
regarding  the  problems  arising  by  the  set  up  of  these  zones  are 
provided  in  Chapter  4.  The  other  two  subroutines  use  a  routine 
(QDAG)  to  do  numerical  integration  only.  Regarding  the  ternary 
diffusion  results,  the  software  described  above  was  applied  only  to 
determine  the  Matano  interface  location  for  individual  elements.  This 
Fortran-77  program  was  run  in  the  VAX  System  at  the  University  of 
Florida. 


CHAPTER  4 

DIFFUSIVITY  EVALUATION  IN  THE  BINARY  Nb-Al  SYSTEM 


This  chapter  presents  the  results  and  their  discussion  of  the 
binary  diffusivity  of  A1  in  b.c.c.  Nb.  The  chapter  is  divided  into  four 
sections.  The  first  one  discusses  the  procedure  performed  for  the 
preparation  of  the  diffusion  couples  used  in  this  part  of  the  study. 
Then  the  results  obtained  for  the  evaluation  of  the  binary  diffusivity 
are  presented.  The  third  section  discusses  these  results,  and  a 
summary  of  the  findings  is  presented  at  the  end  of  the  chapter. 

4.1  Diffusion  Couple  Preparation 

This  study  is  based  on  the  evaluation  of  diffusivity  by  using 
the  traditional  method  of  the  diffusion  couple  whose  fabrication 
began  by  preparing  buttons  of  Nb-4.5  at%  A1  alloy  by  arc-melting  of 
high  purity  materials  (Al,  Nb).  In  order  to  eliminate  macro-  and 
microsegregation  from  the  arc-melted  button,  it  was  deformed 
plastically  by  compression  and  then  annealed  at  1400  °C.  Once 
homogenized,  the  specimen  was  prepared  metallographically  to  be 
bonded  by  the  window-frame  technique  to  a  previously  prepared 
pure  Nb  sample.  This  section  discusses  the  experimental  details 
involved  in  the  preparation  of  the  diffusion  couples. 
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The  dendritic  microstructure  of  the  as  arc-melted  Nb-Al  alloy 
is  shown  in  Figure  4.1.  An  oriented  columnar  grain  structure  permits 
to  assert  the  direction  of  solidification  (from  bottom  to  top).  The 
average  width  of  grains  in  the  center  of  the  samples  was  estimated 
to  be  300  pm .  To  have  such  large  grains  to  start  with  is  helpful  for 
bulk  diffusion  studies  as  it  decreases  the  possible  contribution  of 
grain  boundary  diffusivity. 

The  levels  of  both  macro-  and  micro-segregation  of  the  A1  in 
the  arc-melted  button  were  evaluated  using  micrbprobe  analysis. 

The  macro-segregation  of  aluminum  was  measured  using  a  large 
beam  of  20  pm  in  diameter,  which  was  scanned  over  the  sample  from 
top  to  bottom  using  200  pm  steps.  Figure  4.2  shows  the  resultant 
profile.  A  second  type  of  analysis  was  conducted  to  evaluate  the 
interdendritic  microsegregation  of  aluminum.  This  information  was 
used  to  predict  the  annealing  time  needed  to  eliminate  the 
microsegregation  of  A1  at  the  annealing  temperature.  Figure  4.3 
depicts  an  example  of  the  composition  profile  obtained  from 
dendritic  arms.  In  accordance  with  above  analyses,  the  aluminum 
concentration  in  the  sample  after  the  arc-melting  procedure  was 
generally  around  5.5  atomic  percent  in  the  Al-rich  dendritic  zones. 

The  interdendritic  spacing  was  measured  in  the  middle  of  the 
sample  over  a  structure  similar  to  the  one  shown  in  Figure  4.3a.  The 
average  interdendritic  spacing  was  found  to  be  20  pm  (X=  10  pm).  The 
homogenization  time  at  1400  was  estimated  using  the  data 
regarding  the  interdendritic  spacing,  and  the  A1  concentration 
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Figure  4.1  Optical  micrograph  showing  the  dendritic 

microstructure  in  a  Nb-4.5A1  alloy  after  arc- 
melted.  Note  the  big  columnar  grains. 
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Figure  4,2  A  typical  experimental  composition  profile  for 

aluminum  obtained  by  microprobe  analysis 
showing  the  microsegregation  present  in  an  arc- 
melted  Nb-4.5A1  specimen. 
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Figure  4.3  a)  BSE  micrograph,  and  b)  a  typical 

experimental  composition  profile  for  A1 
obtained  by  microprobe  analysis  showing 
microsegregation  through  a  dendrite  in 
an  arc-melted  Nb-4.5A1  specimen. 
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gradient  in  the  dendrite,  arms  in  conjunction  with  the  A1  diffusion 
coefficient  reported  in  Reference  [6].  Assuming  that  the  concentration 
dependence  of  A1  on  distance  obeys  a  simple  harmonic  behavior  [58], 
the  homogenization  time  for  the  arc-melted  Nb-4.5at%Al  was 
calculated  to  be  8.5  hours  at  1400  ®C. 

b.  Homogenized  Structure 

In  order  to  reduce  the  segregation  observed  in  the  as-cast 
structure,  and  therefore  reduce  the  homogenization  time,  the  arc- 
melted  samples  were  deformed  plastically  in  compression  to  a 
thickness  of  about  2  mm  from  an  original  thickness  of  5  mm.  The 
mechanical  deformation  applied  broke  down  the  dendritic  structure, 
making  it  easier  to  homogenize  the  sample.  In  addition,  the  specimen 

size  and  shape  factors  were  not  a  problem,  since  an  adequate 

specimen  for  bonding  was  easily  obtained  by  trimming  the 

compressed  sample  after  being  annealed. 

It  should  be  mentioned  that  initially,  based  in  the  reported  Nb- 
A1  phase  diagram  [10]  a  composition  of  Nb-7at%Al  was  chosen  for 
binary  diffusion  studies.  This  composition  was  close  to  the  solubility 
limit  of  A1  at  low  temperatures.  After  homogenization  the 
segregation  disappeared.  However,  a  precipitate  had  formed  at  the 
grain  boundaries  (Figure  4.4).  Microprobe  analysis  showed  that  this 
precipitate  has  a  composition  corresponding  to  NbsAl  (17.3  at%), 

whereas  the  matrix  composition  was  6.8  at%.  In  accordance  with  the 
Nb-rich  part  of  the  Nb-Al  phase  diagram  (Figure  2.2),  at  1400  °C  the 
solubility  of  A1  in  Nb  is  about  11%  and  the  second  phase  must  have 
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Figure  4.4  BSE  micrograph  of  the  levitated  and  cooper- 

chilled  quenched  Nb-7.1A1  alloy  after 
homogenization  at  1400®C  for  2  hours  and 
furnace  c?o'ing.  Presence  of  panicles  NbaAl  in 
the  matrix  ((J  phase)  is  observed. 
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precipitated  during  cooling  to  room  temperature  Therefore,  a  lesser 
amount  of  aluminum  was  added  for  the  preparation  of  the  arc¬ 
melting  samples  to  ensure  complete  solubility  of  A1  in  the  |3-Nb 
matrix.  The  final  A1  content  in  alloy  Po  used  was  4.5  at%. 

The  plastically  deformed  samples  with  an  A1  content  of  4.5  at% 
were  homogenized  at  different  temperatures  and  times.  When  the 
temperature  and  time  of  the  heat  treatment  exceeded  1400  °C  and  8 
hours,  respectively,  extensive  precipitation  of  a  second  phase  in  the 
P-Nb  matrix  was  found.  This  new  precipitate  was  observed  not  only 
at  the  grain  boundaries,  but  also  inside  the  grains  of  the  p  phase 
(Figure  4,5),  and  it  had  a  different  morphology  than  that  exhibited 
by  the  Nb3Al.  All  these  characteristics  led  us  to  think  that  the  new 
second  phase  could  be  alumina  AI2O3.  In  order  to  prove  this,  EDX  and 
Auger  analyses  were  conducted  on  a  heat  treated  specimen. 

EDX  spectra  from  the  matrix  and  precipitates,  both  at  grain 
boundaries  and  inside  grains,  were  obtained.  Results  are  presented  in 
Figures  4.6  and  4.7.  The  higher  contents  of  aluminum  and  oxygen  in 
the  precipitates  in  comparison  to  the  matrix  content  suggest  the 
presence  of  aluminum  oxide.  The  larger  Nb  peak  in  the  spectrum  of 
the  precipitate  inside  the  grain  comparing  to  that  at  the  grain 
boundary  is  probably  due  to  the  smaller  thickness  of  the  latter  in 
this  instance. 

The  results  from  the  Auger  analysis  of  these  precipitates 
confirmed  that  they  are  alumina  particles.  Figure  4.8  shows  the 
results  obtained  from  a  particle  which  was  initially  sputtered  to 
remove  the  surface  layer.  The  sputtering  permitted  to  achieve  a 
chemical  analysis  over  a  ’fresh’  surface  approximately  50  A  below 
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Figure  4.5  BSE  micrograph  showing  the  two-phase 

structure  obtained  in  a  Nb-4.5A1  alloy  alter 
annealing  at  1400®C  for  8.5  hrs. 
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c) 


Figure  4.6  (Continued)  c)  EDX  spectrum  obtained  from  a 

precipitate  formed  inside  a  grain  of  a  Nb- 
4.5 AI  alloy  after  annealing  at  1400®C  for  8.5 
hours. 


10 


a) 


10  tim 


Cursor  I  O.06I0K«V  •  0 


t  i 

r 

j 

r 

i  1 

1  1 

I 

L 

1 

i 

( 

; 

i 

i 

■ 

I 

i 

: 

O 

1 

1 

N 

i 

! 

! 

Jl 

iH, 

U 

1 

b)  0.000  vrs  -  8192  S.120 


a)  Secondary  Electron  micrograph,  and  b)  EDX 
spectrum  obtained  from  a  precipitate  formed  at 
the  grain  boundary  of  a  Nb-4.5A1  alloy  after 
annealing  at  1400°C  for  8.5  hours. 
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Figure  4.8  Auger  analysis  of  a  precipitate  located  inside  a 

grain  of  Nb-4.5A1  alloy  annealed  at  1400°C  for 
8.5  hours. 
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the  original  surface,  where  the  particle  is  free  of  contamination  from 
atmospheric  compounds  (i.e.  O2,  CO,  CO2)  that  could  influence  the  final 
result.  When  the  Auger  spectrum  in  Figure  4.8  is  compared  with  that 
of  standard  Auger  alumina  sample  (Figure  4.9)  [73]  an  energy  shift 
of  40  eV  of  both  elements  is  detectable.  The  only  difference  between 
the  two  spectra  is  the  excitation  energy  (3  and  10  keV,  respectively) 
which  only  enhances  the  amplitude  of  the  Auger  signal,  but  does  not 
produce  a  change  in  its  kinetic  energy.  This  displacement  is  due  to 
the  insulating  properties  of  the  alumina  particles,  which  are 
embedded  in  the  P-Nb  matrix.  These  particles  are  charged  by 

electron  beam  during  analysis,  causing  the  kinetic  energies  of  the 
ejected  Auger  electrons  to  increase.  As  for  the  P-Nb  matrix.  Figure 
4.10  shows  the  corresponding  Auger  spectrum.  In  this  case,  no 
displacement  for  any  peak  is  observed.  Therefore,  it  was  concluded 
that  the  resultant  precipitate  after  the  homogenization  treatment 
was  due  to  the  internal  oxidation  of  aluminum  at  the  annealing 
temperature. 

The  homogenization  time  and  temperature  was  chosen  based 
on  a  compromise  between  the  process  of  internal  oxidation  and 

elimination  of  microsegregation.  The  annealing  condition  that 
produced  the  best  result  was  that  involving  a  temperature  of  1400  °C 
for  6  hours.  Figure  4.11  compares  the  microstructures  obtained  by 
annealing  at  1400  °C  for  8.5  and  5  hours.  It  was  also  observed  that 
the  presence  of  oxides  is  more  pronounced  near  the  surface  than  in 
the  middle  of  the  sample  as  shown  in  Figure  4.12.  As  mentioned 

earlier,  the  surface  layers  of  the  sample  were  trimmed  off  before  the 

diffusion  bonding  treatment.  The  presence  of  alumina  at  the  center 
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Figure  4.10  Auger  spectrum  obtained  from  the  matrix  (p  phase)  of  a 
Nb-4.5A1  alloy  after  annealing  at  1400  ®C  for  8.5  hours. 
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a) 


Optical  micrographs  showing  the  microstructure  of 
a  Nb-4.5A1  alloy  after  two  different  homogenization 
treatments:  a)  1400  ®C  for  8.5  hours,  b)  1400  ®C  for 
6  hours. 


Figure  4.11 
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Figure  4.12 


A  low  magnification  BSE  micrograph  of  a  Nb- 
4.5 A1  alloy  after  being  heat  treated  at  1400°C 
for  6  hours.  Note  how  the  amount  of  alumina 
precipitation  lessen  from  top  (sample  surface) 
to  bottom  (sample  middle  section). 
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of  specimens  was  very  scarce,  with  an  average  separation  larger  than 
100  microns  (Figure  4.13).  As  will  be  confirmed  later,  this  separation 
between  alumina  particles  is  much  larger  than  the  penetration  depth 
yielded  by  the  diffusion  treatments.  Also  the  aluminum 
concentration  gradient  in  the  matrix  close  to  the  alumina  particles 
was  very  small,  as  illustrated  by  the  microprobe  analysis  shown  in 
Figure  4.14.  Therefore,  even  though  we  could  not  retard  the  internal 
oxidation  of  the  Nb-Al  alloy,  the  extent  of  it  was  reduced 
considerably  to  such  a  grade  that  was  considered  not  to  have  any 
influence  on  the  diffusivity  measurements.  The  grain  size  achieved 
after  the  homogenization  treatment  was  very  large,  approximately 
ASTM  2  size  (Figure  4.11). 

c.  Diffusion  Bonding 

After  being  homogenized  and  metallographically  prepared,  the 
p -alloy  specimen  as  well  as  the  pure  Nb  sample  were  bonded 
together  in  a  hot-press  furnace.  Initially  chopped  zirconia  (Zr02) 
particles  of  less  than  1  pm  in  size  were  used  as  inner  markers.  Figure 
4.15a  shows  the  resultant  bonded  structure,  which  looks  perfectly 
clean,  free  of  any  kind  of  precipitation.  After  the  diffusion  treatment, 
extensive  precipitation  of  alumina  appeared  near  the  interface  on  the 
Nb-Al  side  as  shown  in  Figure  4.15b.  According  to  thermodynamic 
data  given  in  Reference  [74],  zirconia  and  alumina  have  similar  free 
energy  of  formation  for  equal  moles  of  oxygen.  Alumina  is  more 
stable  at  lower  temperatures  but  Zr02  becomes  more  stable  at  higher 
temperatures.  According  to  the  results  found  in  this  study 
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Figure  4.13  BSE  micrograph  of  a  Nb-4.5A1  sample  as- 

homogenized  condition.  Alumina  presence  in 
the  p-phase  matrix  is  very  scarce. 
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Figure  4.15  A  low  magniHcation  BSE  micrograph  of  a  diffusion 
couple  with  zirconia  markers:  a)  Before  diffusion 
treatment,  and  b)  after  diffusion  treatment  at  1400 
®C  for  4  hours. 


apparently  a  reaction  between  A1  and  Zr02  took  place,  and  therefore 
the  use  of  zirconia  markers  was  discarded  and  yttiria  Y2O3  markers 
were  used  instead.  As  in  the  previous  case,  a  clean,  precipitate -free 
bond  was  obtained  under  the  same  experimental  conditions.  Figure 
4.16  shows  the  structure  at  the  interface  along  with  the  composition 
profile  obtained  using  microprobe  analysis.  The  analytical  result 
suggests  that  the  boundary  is  sharp,  indicating  very  limited 
diffusivity  of  aluminum  in  the  Nb  side,  a  highly  desirable  feature  for 
the  initial  bonding  of  the  diffusion  couples. 

4.2  Diffusivity  Analysis 

The  experimental  concentration  profile  obtained  by  microprobe 
analysis  from  each  of  the  treated  diffusion  couples  were  fitted  to  a 
smooth  curve.  Determination  of  the  Matano  plane  position  permitted 
to  compute  the  interdiffusion  flux,  concentration  gradient  and 
interdiffusion  coefficients  at  each  diffusivity  temperature.  Under  the 
assumption  that  DArNNb»DNb  NAi,  the  intrinsic  diffusivity  of  A1  in  Nb 
as  a  function  of  composition  was  determined.  Then,  the  temperature 
dependency  of  the  coefficient  of  diffusion  of  Al  in  Nb,  as  well  as  the 
corresponding  activation  energy  and  frequency  factor  values  were 
calculated. 
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Figure  4.16 
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4.2.1  Postdiffusion  Microstructural  Analysis  of  Binary  Diffusion 

Couples 

The  use  of  yttria  particles  as  markers  in  the  diffusion  couples 
inhibited  the  alumina  precipitation  along  the  bonded  interface  after 
the  diffusion  heat  treatment.  Figure  4,17  shows  the  resulting 
microstructure  of  a  couple  after  interdiffusion  treatment.  The  alumina 
in  the  structure  of  alloy  Po  was  present  scarcely  not  only  along  the 
bonded  interface  of  the  couples,  but  also  in  the  matrix.  This  fact  helps 
to  explain  that  zirconia  particles  were  the  source  of  oxygen  for 
extensive  formation  of  alumina  near  the  bonded  interface,  but  the 
material  itself  also  had  a  propensity  to  internal  oxidation  of 
aluminum  in  the  alloy.  Chemical  analysis  of  the  interstitial  content  of 
the  alloy  po  after  diffusion  treatment  is  reported  in  Table  4.1, 
together  with  that  for  the  ternary  alloys. 

4.2.2  Measurement  of  Concentration  Profiles 

The  microprobe  analysis  performed  used  a  beam  size  of 
approximately  one  micron.  Therefore,  the  data  for  the  composition 
profile  were  obtained  from  points  one  micron  apart.  A  statistical 
resolution  (standard  deviation)  of  0.1  at.%  in  the  counting  process  was 
obtained  by  employing  a  beam  intensity  of  100  nanoamper  and  a 
dwell  time  of  200  seconds  maximum  at  each  point.  When  either  a 
lower  cunent  or  a  shorter  counting  time  was  employed  the  tail  of  the 
aluminum  concentration  profile  on  the  Nb  side  of  the  couple  was 
missed  (Figure  4.18).  The  compositional  resolution  obtained  by  the 


1  15 


a) 


5D  |im 


Figure  4.17  Backscattered  electron  micrographs  showing 

the  resultant  microstructures  of  a  Nb/Nb-4.5A1 
diffusion  couple  after  diffusion  treatment, 
a)  Near  interface,  b)  Away  from  interface. 
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Table  4.1  Interstitial  content  of  alloys  (*) 

(ppm) 

— 

1 

CO 

1 

O 

310 

420 

1 

2 

130 

370 

62 

U 

58 

1 

48 

Alloy 

O 

CO. 

f'l 

C2. 

c2l 

so 

c 

TS 

U 


u 

s 

>v 

T3 

h* 

>> 

JD 

O 

•o 

:s 

E 


Vl 

R> 

C 

3 


CS 

O 

1 

u 

JZ 

U 


Composition 


1  1  7 


Distance 

(pm ) 

Figure  4.18  A  typical  experimental  composition  profile  for 
aluminum  obtained  by  microprobe  analysis  of  a 
diffusion  couple  Nb/Nb-4.5Al  after  annealing  at 
1450  ®C  for  4.5  hours.  Note  the  sudden  drop  in 
the  A1  composition  profile  in  the  Nb  side  of  the 
couple. 
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microprobe  analysis  is  governed  by  two  factors:  the  background/peak 
intensity  ratio,  and  the  counting  statistics.  When  a  standard  deviation 
larger  than  0.1%  is  employed,  e.  g.  0.25%,  the  counting  lime  for  the 
analysis  is  shorter,  and  the  number  of  events  (x-ray  counts)  detected 
will  be  smaller.  Consequently,  in  regions  where  the  A1  composition  is 
low,  the  x-ray  counts  will  be  insufficient  to  be  taken  into  account,  and 
therefore  the  aluminum  concentration  at  those  points  will  be  detected 
as  zero%.  The  change  to  the  condition  mentioned  above  permitted  to 
detect  the  missing  part  of  the  profile  as  more  counts  were  detected  at 
longer  times.  On  the  other  hand,  it  was  observed  that  employing 
longer  analysis  time  results  in  a  severe  contamination  of  the  surface 
of  the  sample,  as  observed  in  Figure  4.19.  However,  comparison  of  the 
experimental  composition  profiles  obtained  under  the  two  conditions 
revealed  no  differences  in  the  general  result.  The  measurement  of  the 
A1  concentration  at  each  point  took  approximately  6-7  minutes.  In 
average,  to  obtain  an  entire  profile  took  5  hours.  A  minimum  of  two 
profiles  were  measured  for  each  diffusion  couple.  A  typical 
composition  profile  is  shown  in  Figure  4.19. 

4.2.3  Concentration  Profile  Analysis 

The  experimental  raw  data  consisting  of  a  set  of  positions  along 
the  diffusion  direction  and  the  corresponding  concentrations  of 
aluminum  (atomic  %)  were  fitted  to  a  smooth  curve  using  a  modified 
spline  cubic  method  [8].  To  do  so,  generally,  the  set  of  n  data  points  is 
divided  into  intervals  separated  by  knots.  These  knots  set  imaginary 
regions  along  the  composition  profile  in  accordance  to  the  general 
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Figure  4.19  a)  Secondary  electron  micrograph,  and  b)  a 
typical  experimental  composition  profile  for 
A1  obtained  by  microprobe  analysis  of  a 
Nb/Nb-4.5A1  diffusion  couple  after  , annealing 
at  1350  ®C  for  5.5  hours. 
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shape  of  the  curve  in  that  region  (linear,  parabolic,  or  spline  cubic). 
The  spline  cubic  zone  is  also  subdivided  in  intervals  separated  by 
knots  and  then  a  cubic  spline  function  fits  a  different  cubic 
polynomial  to  the  data  in  each  interval  of  this  zone.  For  example,  the 
polynomial  for  the  F*  interval  will  be: 

Cai(x)  =  Cii  +  ci2(x-xi)  +  Ci3(x-xi)2  +  Ci4(x-xi)3  (4.6) 

where  Cai(x)=  Fitted  concentration  curve  in  the  i‘**,  interval 
Cij  =  Cubic  coefficients  for  the  i**'.  interval 
Xi  =  position  of  the  i*  data  point 
Xj  =  position  of  the  i*  knot 

The  coefficients  Cy  and  the  final  knot  positions  are  determined  by 
minimizing  the  least  squared  error  of  the  fitted  curve  to  the 
constraint  that  the  function  Cai(x)  and  its  first  two  derivatives  must 
be  continuous  at  each  knot.  The  number  of  knots  used  for  the  cubic 
spline  fitting  is  important.  The  larger  number  of  knots  used  tends  to 
minimize  the  least  square  error,  i.  e.,  increase  the  goodness  of  the  fit. 
In  addition,  five  data  points  per  interval  must  be  set  so  that  a  third 
degree  polynomial  is  adequately  fitted.  For  the  case  of  the  modified 
spline  cubic  routine  used  in  this  research,  the  domain  of  data  of  the 
diffusion  experiments  are  divided  into  five  different  regions,  as 
shown  in  Figure  4.20.  The  central  part  of  the  data  is  fitted  by  a  cubic 
spline  approximation  as  mentioned  above,  and  the  points  deemed  to 
be  outside  of  the  diffusion  zone  are  fitted  with  horizontal  lines  whose 
ordinates  are  the  end  concentrations  of  the  alloys  forming  the 
diffusion  couple.  The  reason  for  this  horizontal  line  sections  is  that 
cubic  spline  fits  do  not  easily  yield  flat  sections  typical  of  diffusion 
data  outside  the  diffusion  zone.  Since  neither  the  value  of  the  spline 
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Figure  4.20  Schematic  representation  of  the  modified  cubic  spline 
intervals  employed  for  curve  fitting. 
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cubic  nor  its  first  derivative  may  be  specified  at  the  edges  of  the 
region  fitted  with  cubic  splines,  small  parabolic  sections  are  used  to 
smoothly  join  the  cubic  spline  region  with  the  linear  zones.  These 
parabolic  sections  are  constrained  to  match  ordinates  and  slopes  with 
the  horizontal  line  sections  and  ordinates  with  the  cubic  spline  curves. 
When  the  limits  for  the  parabolic  section  were  not  properly  chosen, 
the  fitted  curve  exh  ^ited  small  discontinuities  that  do  not  belong  to 
the  data  (Figure  4.21),  To  avoid  this  problem  several  intervals  for  the 
parabolic  limits  where  tried.  This  matter  is  discussed  in  detail  in  a 
subsequent  section. 

The  program  used  for  curve  fitting  produced  two  output  files. 
One  file  contained  the  polynomial  equations  employed  for  the  other 
subroutines  to  calculate  parameters  such  as  the  Matano  interface, 
interdiffusion  fluxes,  concentration  gradients,  and  coefficients  of 
interdiffusion;  the  other  output  file  contained  100  data  points  of  the 
fitted  curve,  which  were  also  computed  by  the  polynomial  equations 
in  the  previous  file.  Figure  4.22-4.26  contrast  the  experimental 
profiles  with  those  after  curve  fitting  for  the  different  diffusion 
treatments  performed  in  this  investigation.  It  can  be  observed  that 
there  are  not  discontinuities  along  the  entire  profile  and  a  smooth 
curve  could  be  obtained. 

4.2.4  Calculation  of  the  Interdiffusion  Coefficient 

a.  Determination  of  the  Matano  Interface 

The  position  of  the  Matano  interface  was  calculated  in  order  to 
evaluate  the  coefficient  of  interdiffusion  in  the  binary  alloy  couple 
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Figure  4.21  A  typical  composition  profile  for  aluminum 
after  incorrect  spline  cubic  fitting.  Note  the 
discontinuities  along  the  curve. 
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Figure  4.22  Fitting  of  experimental  composition  profiles  by 
the  spline  cubic  method  for  the  Nb/Nb-4.5A1 
couple  heat  treated  at  1350  °C: 

a)  Experimental  curve, 

b)  Fitted  profile. 
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Figure  4.23  Fitting  of  experimental  composition  profiles  by 
the  spline  cubic  method  for  the  Nb/Nb-4.5A1 
couple  heat  treated  at  1400  ®C: 

a)  Experimental  curve, 

b)  Fitted  profile. 
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Figure  4.24  Fitting  of  experimental  composition  profiles  by 
the  spline  cubic  method  for  the  Nb/Nb-4.5A1 
couple  heat  treated  at  1450  ®C: 

a)  Experimental  curve, 

b)  Fitted  profile. 
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Figure  4.25  Fitting  of  experimental  composition  profiles  by 
the  spline  cubic  method  for  the  Nb/Nb-4.5A1 
couples  heat  treated  at  1500  ®C: 

a)  Experimental  curve, 

b)  Fitted  profile. 
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Figure  4.26  Fitting  of  experimental  composition  profiles  by 
the  spline  cubic  method  for  Nb/Nb-4.5A1 
couple  heat  treated  at  1550  °C: 

a)  Experimental  curve, 

b)  Fitted  profile 

Note:  Position  of  the  marker  plane  was  not 
determined  accurately  at  1550  °C;  instead, 
position  of  the  Matano  Plane  in  the  fitted 
profile  is  provided. 
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using  the  Matano-Boltzmann  technique.  The  position  of  the  Matano 
plane  can  be  determined  as  a  function  of  known  quantities  as: 

Xni=(x-CAr  -  x+Cai-^  +  {A})/(Cai-  -  Cai+)  .  (4.13) 

where  {A}  =  entire  area  under  the  curve. 

X-,  x+  =  terminal  position  values  at  the  end  of  the  diffusion  zone. 

Cap*  Cai'*’  =  terminal  composition  values  at  the  end  of  the 
diffusion  zone. 

Equation  4.13  was  applied  to  the  fitted  composition  profiles  computed 
for  A1  concentration  data  from  each  diffusion  couple  as  those  shown 
in  Figures  4.22-4.26. 

b.  Evaluation  of  Jjj^xl 

The  interdiffusion  flux  Ji  was  calculated  as: 

x+ 

Ji(x)  =l/2t  [  (X  -  x^)  Cj(x)  -  (  x+  -  x^)  C  +  +  j  Cjdx  ]  .  (4.16) 

X 

Here  t  is  the  diffusion  time  and  C,+  is  the  terminal  concentration  of 
the  iih  component  at  the  right  end  of  the  diffusion  zone.  Once  the 
Matano  interface  is  evaluated,  all  the  parameters  in  equation  4.16  are 
known,  and  the  interdiffusion  flux  can  be  evaluated  at  each  point  on 
the  concentration  fitted  curve.  Figure  4.27  and  4.28  show  examples  of 
the  computed  interdiffusion  flux  profile,  expected,  the  maximum 
value  of  flux  is  obtained  at  the  Matano  plane  location,  and  it  is  larger 
at  higher  temperatures. 
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a)  Typical  calculated  interdiffusion  flux  profiles 
for  binary  diffusion  couples  Nb/Nb-4.5A1 
interdiffused  at  1350  ®C. 
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Figure  4.28  b)  Typical  calculated  interdiffusion  flux  profiles 

for  binary  diffusion  couples  Nb/Nb-4.5A1 
interdiffused  at  1550  ®C. 
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c.  Evaluation  of  Concentration  Gradient 

The  composition  gradient  was  calculated  by  finding  the 
analytical  derivatives  of  the  fitted  equations  to  the  concentration 
profile.  It  was  found  that  the  calculated  concentration  gradient  was 
very  sensitive  to  the  choice  of  the  limits  for  connecting  the  parabolic 
zone  to  the  linear  and  spline  cubic  regions  during  the  curve  fitting 
procedure.  Even  though  the  concentration  Htted  curve  appeared  to  be 
smooth  and  continuous,  the  fitting  method  did  not  provide  a 
continuous  derivative,  and  discontinuities  and/or  jumps  in  both 
gradient  concentration  and  diffusivity  curves  were  observed  (Figure 
4.29).  It  was  noted  that  the  number  of  knots  employed  into  the  cubic 
spline  region  affected  such  profiles.  By  try  and  error  suitable  limiting 
positions  of  the  different  regions  along  with  different  number  of 
knots,  that  resulted  in  a  smooth  curve  for  the  concentration  gradient 
as  a  function  of  position,  were  obtained.  It  was  noted  that  the  linear- 
parabolic  connection  was  easier  to  fit  than  the  parabolic-spline  cubic 
one.  The  link  of  the  three  regions  was  easier  to  get  in  the  pure  Nb 
side  of  the  couple  as  the  concentration  proHle  was  steeper  at  very  low 
A1  concentrations,  where  the  parabolic  region  extended  for  a  few 
microns.  The  connection  on  the  other  side  of  the  couple,  on  the  other 
hand,  was  much  more  difficult  as  a  larger  penetration  of  the  diffusion 
zone  in  the  Nb-Al  alloy  (Po)  was  obtained,  producing  a  less  step 
concentration  gradient.  This  asymmetry  in  the  concentration  profile, 
with  a  deeper  penetration  on  the  high  aluminum  side,  reflects  a 
difference  in  the  intrinsic  diffusivity  of  the  two  components, 
indicating  that  the  aluminum  may  move  faster  than  niobium. 
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Figure  4.29  Typical  calculated  composition  gradient  and 
interdiffusion  coefficient  profiles  illustrating 
false  peaks  due  to  a  bad  designation  of 
parabolic  regions  during  the  fitting  procedure 
for  binary  diffusion  couples  Nb/Nb-4.5A1 
interdiffused  at  1550  ®C. 

a)  Derivative 

b)  Coefficient  of  Interdiffusion 
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d.  Evaluation  of  Diffusivitv.  Activation  energy,  and  Frequency  Factor 

When  adequate  links  were  set,  smooth  profiles  of  interdiffusion 
coefficient  were  obtained.  Figures  4.30-4.34  present  the  calculated 
composition  gradient  and  composition-dependent  interdiffusion 
coefficient  curves  for  the  b.c.c.  solid  solution  of  the  binary  Nb-Al 
system  in  the  temperature  range  of  1350  ®C  to  1550  °C.  The 
interdiffusion  coefficients  are  summarized  in  Table  4.2,  where  the 
variation  of  the  interdiffusivity  with  the  composition  is  evident  as  it 
drops  at  lower  A1  contents.  The  interdiffusion  coefficient  showed  a 
linear  dependency  on  the  A1  concentration  as  shown  in  Figure  4.35. 
The  diffusivity  coefficient  doubles  its  value  at  a  given  temperature 
with  an  increase  in  concentration  from  1  to  3  at%  Al.  This  strong 
dependency  suggests  that  the  activity  coefficient  of  Al  depends 
considerably  on  the  composition  of  Al. 

The  variation  of  the  interdiffusion  coefficient  as  a  function  of 
temperature  was  found  to  obey  an  Arrhenius  type  behavior.  The 
activation  energy  and  frequency  factor  values  were  computed 
according  to: 

D=  Do  exp  (-Q/rt)  .  (4.17) 

Table  4.3  shows  the  resultant  values  of  Q  and  Do  as  a  function  of 
composition.  The  average  values  for  the  activation  energy  Q  and 
frequency  factor  Do  in  the  temperature  range  from  1350°  to  1550  °C 
were  found  to  be: 


Q=  80.165  Kcal 
Do=0. 1 1 36cm2/sec 
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Figure  4.30  a)  Composition  gradient,  and  b)  interdiffusion 
coefficient  profiles  for  Nb/Nb-4.5A1  diffusion 
couple  heat  treated  at  1350  ®C. 


Intardifluslon  Co«tlici«nt  (E^ll)  ^  Cenc«nlr«tlon  Qr«di«nt 
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Figure  4.31  a)  Composition  gradient,  and  b)  interdiffusion 
coefficient  profiles  for  Nb/Nb-4.5A1  diffusion 
couple  heat  treated  at  1400  ®C. 
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Figure  4.32  a)  Composition  gradient,  and  b)  inierdiffusion 
coefficient  profiles  for  Nb/Nb-4.5A1  diffusion 
couple  heat  treated  at  1450  °C. 
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Figure  4.33  a)  Composition  gradient,  and  b)  interdiffusion 
coefficient  profiles  for  Nb/Nb-4.5A1  diffusion 
couple  heat  treated  at  1500  ®C. 
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Figure  4.34  a)  Composition  gradient,  and  b)  interdiffusipn 
coefficient  profiles  for  Nb/Nb-4.5A1  diffusion 
couple  heat  treated  at  1550  ®C. 
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4.2.5  Calculation  of  Intrinsic  Diffusivities 

In  accordance  with  Darken’s  analysis  (section  2.2. 1.4),  equations 
2.38  and  2.39  permit  to  evaluate  the  intrinsic  diffusion  coefficients  of 
binary  system  when  the  experimental  parameters  such  as  marker 
velocity  and  concentration  gradient  are  known.  The  results  of  this 
research,  however,  indicated  that  the  position  of  the  Matano  interface 
is  very  close  to  the  markers  position;  the  distance  between  the  final 
marker  position  and  the  Matano  interface  for  the  largest  depth  of 
penetration  (1550  ®C,  5  hours)  was  less  than  a  few  micrometers, 
which  is  within  the  accuracy  of  the  measurements.  Therefore,  the 
interface  velocity  could  not  be  calculated  accurately.  Tlie  intrinsic 
diffusivity  of  A1  was  estimated  (Table  4.4)  by  assuming  that 
Dai  NNb  »  DNb  Nai; 

D  »  Dai  NNb 

This  assumptions  is  justified  by  the  fact  that  the  self-diffusion 
coefficient  of  Nb  in  the  temperature  range  of  this  study  is  in  the  order 
of  10"  cm2/sec.  [22,  23].  It  is  acknowledge  that  alloying  increases 
the  diffusivity  of  Nb,  however,  it  is  expected  to  be  orders  of 
magnitude  smaller  than  -10*1 1  cm^/sec,  which  is  the  interdiffusion 
coefficient  calculated  for  the  Nb-Al  binary  system.  Furthermore 
NAl<<NNb  in  the  couples  studied.  Figures  4.36  to  4.40  present  the 
Arrhenius  plots  for  different  aluminum  concentrations  (1%  to  3%). 
The  calculated  activation  energy  and  frequency  factor  values  are 
given  in  Table  4.5.  The  diffusivity  data  for  lower  and  higher  A1 
contents  were  not  considered  since  large  errors  are  obtained  in  the 


Figure  4.36  The  temperature  dependency  of  the  intrinsic 

diffusion  coefficient  of  A1  in  Nb  at  the  aluminum 
concentration  of  3.0  at%. 
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Figure  4.37  The  temperature  dependency  of  the  intrinsic 

diffusion  coefficient  of  A1  in  Nb  at  the  aluminum 
concentration  of  2.5  at%. 
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Figure  4.38  The  temperature  dependency  of  the  intrinsic 

diffusion  coefficient  of  A1  in  Nb  at  the  aluminum 
concentration  of  2.0  at%. 
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Figure  4.40  The  temperature  dependency  of  the  intrinsic 

diffusion  coefficient  of  A1  in  Nb  at  the  aluminum 
concentration  of  1.0  at%. 
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integration  step  when  the  composition  is  close  to  the  terminal 
compositions.  The  linear  curve  fitting  obtained  is  very  good  as  shown 
by  the  values  of  the  correlation  coefficient  R2  (>0.95)  in  the  worst 
case).  The  average  values  of  the  activation  energy  and  frequency 
factor  for  the  diffusivity  of  A1  in  Nb  in  the  composition  range  of  1.0  to 
3.0  atomic  percent  of  A1  are  calculated  to  be 

Do  =  0.1161  cm2/sec, 

Q  =  80.2  Kcal/mole 

4.3  Discussion 

There  has  been  only  one  reported  study  on  the  diffusivity  of  A1 
in  b.c.c.  Nb  [6].  As  detailed  in  Chapter  2,  Bodrov  and  Nikolaev  [6] 
calculated  the  diffusion  coefficient  of  A1  in  Nb  at  the  temperature 
range  of  1500®  to  1950  ®C  by  an  analysis  of  the  kinetics  of 
evaporation  of  A1  present  as  a  dilute  impurity  in  Nb  in  a 
concentration  range  of  10*2- 10*3  at%  using  atomic  absorption 
spectroscopy.  In  order  to  compare  the  results  of  this  study  with 
theirs,  the  Dai  values  were  plotted  as  a  function  of  A1  concentration. 
The  Dai  values  for  dilute  solutions  were  estimated  by  the 
extrapolation  of  the  A1  intrinsic  diffusivity  to  zero  concentration.  With 
the  extrapolated  results  a  new  Arrhe.iius  curve  was  obtained  which  is 
shown  in  Figure  4.41  together  with  their  results.  The  Dai  values 
calculated  in  this  study  are  an  order  of  magnitude  smaller. 
Furthermore  the  activation  energy  reported  by  Bodrov  and  Nikolaev 
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Figure  4.41  Comparison  between  this  work  and  Ref.  [6]  for  the 

temperature  dependence  of  the  diffusion  coefficient 
of  AI  in  Nb. 
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[6]  is  much  larger  (102  Kcal/mole)  than  the  value  obtained  in  this 
study  (74.3  kcal/mole).  These  discrepancies  can  be  discussed  in  terms 
of  the  differences  between  the  two  studies; 

(i)  Method:  The  indirect  method  of  evaporation  of  A1  employed 
by  Bodrov  and  Nikolaev  involves  assumptions,  which  their  validity  is 
not  justified  (see  section  2.1.2).  The  traditional  method  employed  in 
this  sudy  is  a  direct  method  and  is  believed  to  yield  more  accurate 
results. 

(ii)  Temperature  range:  This  study  was  conducted  a  lower 
temperature  range  1350-1550  ®C  in  comparison  to  1500-2000®C, 
which  was  applied  to  their  experiments.  The  self-diffusivity  of  Nb  is 
reported  to  be  very  high  close  to  its  melting  point  [23].  This 
phenomenon  has  been  attributed  to  the  rapid  motion  of  divacancies 
at  high  temperatures  [52].  This  section  was  explained  in  detail  in 
section  2.1.2.  The  contrioution  of  divacancies  at  high  temperatures 
can  explain  a  larger  activation  energy  and  higher  diffusion  coefficient 
at  temperatures  close  to  the  solidus  temperature  of  the  alloy. 
However,  since  there  is  an  overlap  of  test  temperature  between  this 
study  and  that  of  Bodrov  and  Nikolaev  the  low  diffusivity  of  A1 
measured  in  this  study  can  not  be  explained  solely  by  a  divacancy 
mechanism. 

(iii)  Interstitial  Content:  The  processes  of  bonding  and  diffusion 
treatment  introduced  a  high  concentration  of  nitrogen  and  oxygen  in 
the  alloys  used  in  this  study  (see  Table  4.1).  Eizinger  et.  al.  [23]  have 
reported  that  oxygen  content  has  no  effect  on  the  self-diffusion  of  Nb. 
Al  is  expected  to  bind  with  oxygen,  however,  such  trapping  would 
expect  to  increase  the  activation  energy  in  comparison  to  a  purer 
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system.  Therefore,  the  observed  inconsistency  can  not  be  attributed 

to  the  relatively  high  concentration  of  interstitials  in  the  alloys  used 
in  this  study. 

Bodrov  and  Nikolaev  also  measured  the  intrinsic  diffusivity  of 

Al  in  Ta: 

Do=  1.525  cm2/sec 

Q=  73.14  kcal/mol 

Di5oo«=1.3646  X  lO-^cm^/sec 

It  is  interesting  to  note  that  the  Q  value  for  Ta-Al  is  much  smaller  and 
closer  to  the  Q  value  found  in  this  study  for  the  Nb-Al  system. 
Furthermore  a  recent  work  by  Maslov  et.  al.  [75],  who  have 
investigated  the  diffusion  of  Al  in  V  in  the  temperature  range  of 
1000-1500  ®C,  gives: 

Do=  0.084  cm2/sec 
Q=  64.1  kcal/mole 

The  activation  energies  for  self-diffusion  of  Ta,  Nb  and  V  normalized 
for  melting  point  (Q/Tm)  are  very  close  [28,  76]  (see  Figure  2.10),  and 
the  phase  diagram  of  these  elements  with  Al  are  very  similar. 
Therefore  one  would  expect  that  the  activation  energy  for  diffusivity 
of  Al  in  Nb  to  be  in  between  that  in  Ta  and  V,  consistent  with  the 
result  obtained  in  this  study. 

Most  of  the  existing  interdiffusion  studies  including  Nb  as  one  of 
the  components  are  those  involving  another  transition  metal  as  the 
second  component  in  the  diffusion  couple.  Table  4.6  [77]  present  the 
value  of  the  activation  energy  for  the  binary  system  Nb-X  (X=  Ti,  V, 
Zr,  Mo,  Ta,  W,  Al).  In  this  Table  the  intrinsic  diffusivities  of  the  solutes 
are  compared  to  the  self-diffusion  of  Nb.  All  these  elements  form 
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extensive  b.c.c.  solid  solution  regions  in  the  phase  diagram  when 
combined  with  Fb.  Except  for  Mo  and  W,  the  other  elements  show 
relatively  a  similar  activation  energy  for  diffusion.  Although  A1  is  not 
a  transition  metal,  the  results  of  the  ternary  diffusivity  obtained  in 
this  study  (see  Chapter  5)  suggest  that  in  the  Nb  comer  of  the  ternary 
phase  diagram  A1  and  Ti  have  similar  penetration  tendencies  at  1200 
®C  and  1400  “C  annealing  temperatures.  Again,  this  observation 
suggests  that  in  comparison  to  the  Bodrov  and  Nikolaev  study  the 
lower  diffusivity  coefficients  obtained  in  this  study  correlate  better 
with  the  diffusivity  of  Ti  in  Nb. 

The  results  of  interdiffusivity  for  different  phases  in  the  Nb-Al 
system  are  summarized  in  Table  4.7  [78-80].  The  activation  energy 
for  the  interdiffusion  in  the  p  phase  is  somewhat  smaller  than  NbaAl. 
The  activation  energy  seems  to  decrease  with  increasing  the  A1 
content  and  for  NbAla  it  reaches  the  activation  energy  for  self¬ 
diffusion  in  pure  A1  (35  kcal/mole)  [77]. 

4.4  Summary 

The  binary  diffusion  in  the  P  Nb-Al  solid  solution  was 
evaluated  by  the  Matano-Boltzmann  method.  The  Nb-4.5A1  alloy 
used  for  this  research  was  susceptible  to  internal  oxidation.  The 
formation  of  AI2O3  near  the  interface  could  be  lessened  substantially 
by  using  yttria  markers  instead  of  zirconia  ones.  Due  to  the  oxidation 
restriction  both  temperature  and  time  of  interdiffusion  treatment 
had  to  be  shortened  and  therefore,  the  extent  of  composition  profiles 
obtained  was  limited  to  a  maximum  of  45  pm.  As  a  result. 
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appreciable  movement  of  markers  after  interdiffusivity  was  not 
detected  either. 

Spline  cubic  fitting  of  the  experimental  composition  profiles 
worked  successfully,  even  though  difficulties  arose  during  curve 
fitting,  specially  in  the  parabolic  region  of  the  profile  in  the  Po  alloy 
side  of  the  diffusion  couple. 

The  results  of  this  study  indicate  that  the  interdiffusion 
coefficient  in  |5  solid  solution  of  Nb-Al  system  varies  linearly  with  A1 
composition  in  the  concentration  range  of  1-3  at%  A1  and  it  follows  a 
linear  Arrhenius  behavior  as  function  of  temperature  in  the  1350  ®C- 
1500  ®C  range.  The  activation  energy  for  diffusivity  of  A1  in  Nb  was 
found  to  be  80.2  kcal/mole,  wich  is  smaller  than  the  previously 
reported  value  of  102  kcal/mole  that  has  been  obtained  in  the  1500 
°C-1900  ®C  temperature  range.  It  is  suggested  that  this  discrepancy  is 
associated  with  the  differences  in  the  evaluation  method  and 
possibly  the  temperature  range. 


CHAPTER  5 

DIFFUSIVITY  ANALYSIS 
IN  THE  TERNARY  Nb-Ti-Al  SYSTEM 


The  purpose  of  this  part  of  the  study  was  to  investigate  the 
diffusivity  in  B.C.C.  solid  solution  of  ternary  Nb-Ti  A1  system.  The 
ternary  diffusion  couples  prepared  in  this  study  are  illustrated  in 
Figure  5.1. 

The  Nb-42.5at%Ti-15at%Al  (Pi)  alloy  was  chosen  as  the  main 
alloy  and  it  was  coupled  with  pure  Ti,  pure  Nb  and  the  Nb-50at%Ti 
(P)  alloy.  The  composition  of  A1  was  chosen  based  on  the  information 
available  on  the  p-transus  temperature.  An  increase  in  the  A1  content 
in  this  part  of  the  ternary  phase  diagram  results  in  an  increase  in  the 
P-transus  temperature  and  order-disorder  temp'^.ratures  [81]. 
Therefore  15  at%  was  estimated  to  be  the  maximum  A1  content  of  an 

alloy  that  would  have  a  p-transus  and  order-disorder  temperatures 

below  1000  ®C.  Two  additional  couples  consisting  of  the  Nb- 

47.5at%Ti-5at%Al  (P2)  alloy  coupled  with  the  P3  and  Po  alloys  were 
also  prepared  to  complement  the  binary  diffusivity  results. 

In  this  chapter  the  results  of  the  diffusivity  experiments 
performed  on  the  ternary  Nb-Al-Ti  system  and  their  discussion  are 
presented.  The  first  part  of  the  chapter  deals  with  the 

microstructural  details  of  the  ternary  couples  before  and  after  the 
diffusion  treatment.  Then,  in  order  to  have  a  better  understanding  of 
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Figure  5.1  Schematic  rep'^'^sentaiion  of  the  ternary  diffusion  couples  employed  in  this 
research. 
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the  diffusivity  behavior  of  the  system,  the  composition  profiles 
obtained  by  microprobe  analysis  are  given  in  different  groups 
depending  on  the  compositional  characteristic  of  the  couple.  The 
information  provided  by  these  composition  profiles  is  presented  and 
analyzed  with  the  help  of  concepts  such  as  Penetration  Tendencies, 

Composition  Paths,  Kirkendall  Shifts,  and  Crossover  Displacements.  A 
summary  of  the  findings  is  given  at  the  end  of  the  chapter. 

5.1  Microstructural  Analysis 

The  microstructure  of  the  individual  alloys  and  the  diffusion 
couples  after  both  the  bonding  and  the  diffusion  treatments  were 

investigated  using  optical  and  electron  microscopy  techniques. 

Pure  Ti  undergoes  an  allotropical  transformation  at  882  ®C, 
from  a  high  temperature  body-centered  cubic  (P)  phase  to  a 
hexagonal  close-packed  (a)  structure  (see  Figures  2.1  and  2.3).  The  Ti 
alloy  phase  diagrams  in  the  Ti-rich  section  are  usually  classified  into 
three  types.  Elements  that  stabilize  a -phase  and  increase  the  a/p 
transus.  Aluminum  and  oxygen  are  the  most  important  a  stabilizers. 
Elements  that  depress  the  a/p  transus  and  stabilize  p -phase  are 
classified  in  two  groups:  those  which  form  binary  systems  of  the  P- 
isomorphous  type,  such  as  Nb  (see  Figure  2.1)  and  those  which 

favour  formation  of  a  P-eutectoid  such  as  hydrogen.  In  general  the 
microstructure  observed  is  sensitive  to  the  cooling  rate  from  the  P- 
phase  regime.  A  fast  cooling  rate  can  result  in  formation  of 

martensitic  a’  phase. 
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The  as-received  Ti  used  in  this  study  had  a  duplex 
microstructure  consisting  of  a  |3  phases  as  shown  in  Figure  5.2.  The 
small  amounts  of  impurities  (see  Table  3.1)  are  sufficient  to  result  in 
an  (a  P)  microstructure  rather  than  a  single  a -phase  at  room 
temperature.  On  the  other  hand  Pi(Nb-42.5Ti-15Al),  P2  (Nb-47.5Ti- 
5A1)  and  P3  (Nb-50Ti)  alloys  displayed  a  single  phase  microstructure 
as  shown  in  Figure  5.3,  while  the  phase  diagram  suggests  that  these 
alloys  should  consist  of  a  and  P  phases  at  room  temperature.  This 
observation  can  be  explained  by  the  low  mobility  of  Nb  at  the  a/p 
transition  temperature  for  these  alloys.  TEM  analysis  of  a  50Nb- 
40Ti-10Al  alloy  [82]  has  revealed  that  this  alloy  is  composed  of  the 
disorder  P -phase.  Although  the  exact  P-transus  temperature  for 
these  alloys  are  not  known,  based  on  the  binary  Nb-Ti  phase 
diagram,  they  can  be  estimated  to  be  below  1000  ®C. 

As  shown  in  Figure  5.4  the  duplex  microstructure  of  the  pure 
Ti  persisted  after  bonding  with  the  pi  alloy  at  a  bonding  temperature 
of  700  °C,  which  is  below  the  P-transus  temperature  of  titanium. 
However,  in  the  vicinity  of  the  interface,  where  some  interdiffusion 
had  occurred,  the  microstructure  consists  of  a  single  phase  as  shown 
in  Figure  5.5a.  The  microprobe  analysis  showed  slight  penetration  of 
both  A1  and  Nb  to  the  single  phase  region  as  given  in  Figure  5.5b. 
Since  A1  is  a  strong  a-stabilizer  it  may  be  speculated  that  this  narrow 
zone  of  penetration  consists  of  a-phase.  Figure  5.5  also  shows  that  a 
thin  darker  band  has  developed  in  the  ternary  side  of  the  couple. 
This  is  probably  due  to  the  slight  interdiffusion  of  Ti  into  the  P 1 
phase  which  is  detected  by  the  composition  dependent  BSE  image.  A 


b) 


Micrographs  showing  ihe  duplex  structure 
(a+P)  of  the  as-received  Ti:  a)  Optical 
micrograph,  as-etched  condition;  b)  BSE 
image. 


Figure  5.2 
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Ij  Figure  5.3  BSE  micrograph  of  pi  alloy  (Nb-42.5Ti-15Al)  after 

homogenization  at  1400  ®C  for  6  hours. 
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Figure  5.4  Optical  micrograph  of  the  bonded  interlace 

in  diffusion  couple  #12  (pl/Ti)  after  bonding 
at  700  °C  for  2  hours  showing  the  duplex 
structure  a+P  of  the  Ti  side  of  the 
couple. 
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Figure  5.5  a)  BSE  micrograph  of  interface  in  diffusion 

couple  #12  (Pl/Ti)  before  diffusion  treatment; 
b)  the  three  composition  profiles  across  the 
boundary 
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diffusion  treatment  at  1200  ®C  resulted  in  a  two-phase  region 
adjacent  to  the  interface  (Figure  5.6).  The  microstucture  resembles 
the  typical  Widmanstatten  (a  +  P)  structure  that  developes  in  Ti 
alloys.  The  development  of  a  two-phase  region  adjacent  to  the 
interface  in  the  1200  °C  treated  couple  can  be  attributed  to  the  facts 
that  (i)  at  1200  °C  the  microstructure  consists  of  single  p  phase,  (ii) 
the  higher  temperature  allows  the  diffusion  of  Nb  into  the  Ti  side, 
and  (iii)  the  cooling  rate  achieved  in  the  vacuum  furnace  is  faster 
than  that  obtained  in  the  hot  press,  which  was  used  for  the  bonding 
treatment.  The  acicular  morphology  of  this  two-phase  region  in 
comparison  to  the  equiaxed  morphology  of  Ti  can  be  explained  by 
the  higher  alloying  content  of  it. 

Similarly  the  other  ternary  couples  that  were  bonded  at  900 
showed  the  formation  of  alpha  phase,  which  precipitated  along  grain 
boundaries  and  at  the  bonded  interface  as  shown  in  Figure  5.7.  These 
areas  have  a  higher  surface  energy,  which  induced  the  formation  of  a 
phase.  Figure  5.8  shows  BSE  micrographs  where  this  precipitate  looks 
dark  indicating  its  high  A1  content.  The  effect  of  the  aluminum 
content  on  the  presence  of  this  second  phase  is  corroborated  by  the 
following  observation.  No  precipitation  was  observed  in  the  binary 
Ti-Nb  side  of  any  couple  (#7,  #9).  However,  when  a  small  amount  of 
A1  (0.5  at%)  was  present  in  the  Nb-Ti  side  (ps)  of  couple  #  10, 
precipitation  was  observed  on  both  sides  of  the  couple  as  shown  in 
Figure  5.9.  After  interdiffusion  treatment  at  either  1200°  or  1400  °C, 
no  a  phase  was  detected  anywhere  (Figure  5.10),  which  means  that 
due  to  the  higher  temperature  the  structure  transformed  completely 
to  P  phase  and  because  of  the  redistribution  of  A1  and  Nb  around  the 
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Figure  5.6  Optical  micrograph  of  the  Ti  side  next  to  the 

bonded  interface  in  diffusion  couple  #12 
(Pl/Ti)  after  diffusion  treatment  at  1200  °C 
for  6  hours.  Note  the  duplex  (a+p)  zone  with 
the  acicular  morphology  next  to  the  interface 


b) 


BSE  micrographs  showing  precipitation  of 
a  phase  in  the  ternary  side  of  the  couples 
along  the  interface:  a)  diffusion  couple  #9 
(pl/P3);  and  b)  diffusion  couple  #11  (|52/P0). 


Figure  5.8 
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BSE  micrograph  of  bonded  interface  in 
diffusion  couple  #10  (|32/p3)  after  bonding  at 
900  °C  for  2  hours  showing  the  presence  of 
second  phase  (a)on  both  sides  of  the  interface. 


Figure  5.9 
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Figure  5.10  BSE  micrograph  of  the  interface  in  diffusion 
couple  #6  (Pl/Nb)  after  diffusion  treatment 
at  1400  ®C  for  6  hours.  Microprobe  analysis 
contamination  over  the  precipitation-free 
structure  is  detailed. 
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interface  and  the  relatively  faster  cooling  rate,  in  the  vacuum 
furnace  the  precipitation  of  a  phase  from  P  was  inhibited.  After  the 
diffusion  treatment,  the  couples  were  reheated  at  920  °C  in  the 
vacuum  furnace  in  order  to  study  the  stability  of  the  microstructure. 
Again  precipitation  occurred  at  the  grain  boundaries,  interfaces,  and 
on  the  dislocations  substructure,  but  there  was  a  denuded  zone  near 
the  interface  (Figure  5.11),  which  matched  the  diffusion  zone  after 
the  heal  treatment  at  1200  “C.  These  observations,  in  addition  to  the 
fact  that  the  structure  of  ternary  alloys  after  homogenization  at  1400 
consists  of  the  single  f)  phase,  lead  to  the  conclusion  that  at  the 
interdiffusion  temperature  of  1200  ®C  and  above  the  microstructure 
was  composed  of  single  b.c.c.  P  solid  solution. 


5.2  Composition  Profiles 

The  composition  profiles  for  the  different  ternary  diffusion 
couples  after  interdiffusion  at  1200°  and  1400  °C  for  6  hours,  as 
obtained  by  microprobe  analysis  are  shown  in  Figures  5.12  to  5.18. 
Note  that  the  profiles  are  presented  in  different  scales.  At  least  two 
composition  profiles  were  obtained  for  each  couple,  as  it  has  been 
discussed  in  Appendix  A.  Next,  the  composition  profiles  for  two 
different  cases  are  discussed,  and  the  presence  of  critical  points  and 
the  effect  of  temperature  are  presented. 
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Figure  5.11  BSE  micrograph  of  the  p2  side  away  of  the 

interface  of  diffusion  couple  #10  (j32/(53)  after 
diffusion  treatment  at  1200  °C  for  6  hours, 
and  subsequent  treatment  at  920  °C  for  2 
hours.  Reprecipitation  of  the  second  phase  is 
observed. 
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Penetration  curves  of  the  3  components  in  the 
couple  #8  (Pl/Nb)  heat  treated  at  1200  “C  for  6 
hours. 
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The  penetration  curves  of  the  three  couples  containing  Pi  alloy 
are  compared  for  individual  elements  in  Figures  5.19  to  5.21.  It  is 
observed  that  the  depth  of  penetration  in  the  counterpart  side  of  pi 
is  in  general  larger  as  its  Ti  content  increases.  For  example, 
penetration  curve  for  couple  #8  shows  only  a  shallow  penetration  of 
the  three  elements  (less  than  5  jim)  into  the  Nb  side,  whereas  couples 
#9  and  #12  show  much  deeper  penetration  into  the  P3  and  Ti, 
respectively.  As  to  the  penetration  in  the  ternary  alloy  side  Pi,  the 
depth  of  penetration  for  each  element  is  also  larger  as  the  Ti  content 
of  the  counterpart  side  increases. 

5.2.2  Effect  of  Concentration  Gradient 


Figures  5.22  to  5.24  compare  the  composition  profiles  of  each 
element  for  two  cases  with  a  constant  Nb:Ti  ratio  and  different  A1 
composition  gradient  (couples  #9  and  #10).  As  expected,  larger 
composition  gradient  results  in  a  larger  penetration  depth  of  Al. 
When  the  Al  concentration  on  both  sides  of  the  couple  is  the  same 
the  composition  profile  is  given  in  Figure  5.25.  A  larger  diffusion 
zone  on  the  ternary  side  of  the  couple  than  that  of  the  binary  Nb- 
4.5A1  side  is  observed.  However,  when  the  concentration  profiles  for 
diffusion  couples  #11  (P0/P2)  and  #10  (P3/P2)  are  compared  to  each 
other  some  interesting  observations  are  made  (see  Figures  5.26- 
5.28).  It  is  observed  that  although  the  concentration  gradients  of  Ti 
and  Nb  are  smaller  in  couple  #10,  it  has  a  larger  diffusion  zone  than 


2 


204 


the  couple  #1 1  does.  This  could  be  explained  readily  by  comparing 
the  solidus  temperatures  of  the  alloys  involved  in  the  respective 
couples;  the  Po  alloy  has  a  much  higher  solidus  temperature  than  the 
P3  alloy  does. 

5.2.3  Critical  Points 

Except  in  a  few  cases,  all  composition  profiles  of  the  couples 
studied  did  not  exhibit  critical  points  and  their  shape  was  smooth 
and  similar  to  that  obtained  in  the  binary  diffusion  studies.  Figures 
5.29  to  5.31  shows  the  composition  profiles  for  those  cases,  in  which 
maxima  and  minima  where  observed.  An  interesting  observation  is 
that  these  critical  points  appeared  always  when  the  concentration 
gradient  of  the  element  was  small,  i.  e.  for  AI  in  couple  #11  (Figure 
5.29),  for  Nb  in  couples  #7,  #9,  and  #10  (Figure  5.30),  and  for  Ti  in 
couple  #9  (5.31).  When  the  diffusion  temperature  is  increased  for  the 
couple  P3/P1  from  1200  ®C  (#9)  to  1400  °C  (#7)  the  maxima  in  the 
composition  profiles  of  the  Ti  is  completely  removed,  and  for  the  Nb 
at  1400  ®C  the  profile  shows  a  small  tendency  for  maximum.  These 
observations  suggest  that  limited  physical  movement  of  atoms  due  to 
a  low  interdiffusion  temperature,  low  concentration  gradient  or  high 
melting  point  of  the  element  such  as  Nb  results  in  the  occurrence  of 
maxima  and  minima  in  the  concentration  curves  as  obtained  by 
microprobe  analysis.  As  suggested  by  DeHoff  et.  al.  [69]  the  apparent 
diffusion  zone  measured  by  the  microprobe  analysis  does  not 
represent  the  actual  physical  distance  over  which  the  atoms  have 
moved.  When  the  movement  of  one  element  is  limited  due  to  the 
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Figure  5,30  Critical  point  observed  in  the  Nb  composition 
profile  after  diffusion  treatment  of  couples  #7 
(P0/p2),  #9  (p3/pi)  and  #10  (P3/P2)  at  1200  °C 
for  6  hours. 
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Critical  point  observed  in  the  Ti  composition 
profile  after  diffusion  treatment  of  couple  #9 
P3/|31  at  1200  ®C  for  6  hours. 
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above  mentioned  reasons,  such  as  Nb  in  the  Nb-Ti-Al  system  the 
atomic  fraction  profile  of  this  element  represent  the  unequal 
mobility  of  the  other  constituent  elements. 

5.2.4  Temperature  Effect 

The  composition  profiles  for  each  element  of  the  systems 
consisting  of  Nb/Pi/Nb  (#6,  #8),  and  P3/P1/P3  (#7,  #9)  interdiffused 
at  1400  "C  and  1200  ®C  are  compared  in  Figures  5.32  to  5.34  for  the 
former  case  and  Figure  5.35  to  5.37  for  the  latter  one.  At  1400  °C  the 
general  behavior  observed  in  the  diffusion  couples  treated  at  1200  ®C 
is  maintained.  That  is,  the  penetration  depth  is  larger  in  the  ternary 
side  of  the  couples.  As  expected,  the  depth  of  penetration  of  each 
component  becomes  deeper  as  the  temperature  is  raised. 

5.3  Penetration  Depth 

As  was  shown  in  Figures  5.12  to  5.18  the  depth  of  the 
penetration  of  the  composition  profiles  is  not  symmetrical  in  both 
sides  of  a  given  couple.  Table  5.1  shows  the  depth  of  penetration  for 
the  different  couples  as  measured  by  microprobe  analysis.  Next,  each 
case  is  analyzed  individually. 

a)  PjL/Nb  (1200  ®C):In  this  case  the  three  components  show  an  equal 
penetration  depth  on  each  side  of  the  couple;  however,  their 
penetration  curves  are  not  symmetrical.  A  much  larger  penetration  is 
observed  on  the  |3i  side,  i.  e.  the  side  with  the  lower  solidus 
temperature. 
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Figure  5.34 


Composition  profile  for  Ti  after  interdiffusion 
treatment  at  1200  ®C  (#8)  and  1400  "C  (#6) 
for  6  hours. 
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Figure  5.36  Composition  profile  for  Nb  after  interdiffusion 
treatment  at  1200  ‘’C  (#7  )  and  1400  °C  (#9) 
for  6  hours. 
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5.37  Composition  profile  for  Ti  after  interdiffusion 

treatment  at  1200  °C  (#7)  and  1400  °C  (#9)  for  6 


A 


Table  5.1  Depth  of  penetration  in  ternary  diffusion  couples  of  system 

(B.C.C.)  Nb-Ti-Al 
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b)  fi-iZfii  (1200  °C):ln  this  case  the  solidus  temperatures  of  the  two 
sides  are  not  expected  to  be  very  different,  and  therefore  the 
concentration  gradients  of  the  components  play  an  important  role  in 
the  final  shape  of  the  penetration  curve.  Aluminum  is  the  element 
with  the  largest  composition  gradient  and  penetrates  from  the  Pi  to 
the  Ps  side  producing  an  approximately  similar  penetration  depth  on 
both  sides  of  the  couple.  The  flux  opposite  to  the  A1  motion  is  mainly 
due  to  Ti,  as  indicated  by  the  small  penetration  of  Nb  in  the  Pi  side. 

It  is  interesting  that  A1  penetration  depth  in  the  P3  side  is  much 
smaller  than  the  draw  depth  of  Nb  and  Ti.  A  comparison  of  the 
diffusion  zone  size  with  that  at  1400  ®C  for  the  same  couple  reveals 
that  the  apparent  draw  depth  of  Nb  on  the  P3  side  decreases  with 
increasing  temperature.  This  observation  suggests  that  indeed  the 
apparent  draw  depth  of  Nb  at  1200  ®C  does  not  represent  the  actual 
physical  motion  of  Nb.  However,  since  A1  does  not  penetrate  to  the 
same  depth,  unless  the  crystal  can  tolerate  very  large  fractions  of 
vacancies  there  should  be  an  actual  uphill  diffusion  of  Nb.  Another 
possibility  that  can  contribute  to  this  observation  is  a  large 
composition  dependency  of  the  molar  value  in  this  system. 
Unfortunately  no  precise  lattice  paramenter  data  for  these 
compositions  are  known.  Here,  the  Ti  comes  from  a  large  distance 
within  the  (33  side.  The  apparent  large  penetration  depth  of  Nb  on  the 
P3  side  is  due  to  the  readjustment  of  the  atomic  fraction  due  to  the 
lost  of  Ti.  It  should  be  noted  that  the  penetration  depths  evaluated 
by  microprobe  analysis  may  not  represent  the  physical  distance  that 
atoms  have  moved,  since  a  variation  in  their  atomic  fraction  may  be 
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a  result  of  the  diffusivity  of  any  of  the  other  species  in  that  part  of 
the  couple  [69], 

c)  fiLiZIi  (1200®C):  The  relative  penetrations  of  A1  and  Ti  are  equal  on 
each  side  of  this  couple  and  larger  than  that  of  Nb.  It  is  observed 
again  that  the  side  with  the  lower  solidus  temperature  (Ti-side) 
exhibits  the  larger  penetration  depths.  The  depths  of  penetration  in 
this  case  are  in  general  much  larger  as  compared  to  those  in  the 
Pl/Nb  couple.  In  the  latter  the  kinetics  are  limited  by  the  motion  of 
atoms  and  vacancies  in  the  Nb  side,  however,  in  the  former  couple  it 
is  the  motion  in  the  (3i  side  that  controls  the  penetration. 

d)  (1200  ®C):  In  this  case  the  penetration  depth  of  A1  is  smaller 
that  the  Pi/Pa  couple  due  to  a  smaller  A1  concentration  gradient.  It  is 
interesting  to  note  that  in  this  couple,  contrary  to  the  P1/P3  couple 
the  Nb  and  Ti  draw  depths  are  much  smaller  and  equal  to  the 
penetration  depth  of  the  A1  on  the  P3  side.  This  behavior  may  be 
related  to  the  lower  concentration  gradient  of  Ti  in  the  |32/P3  couple 
than  in  the  Pi/^3  couple. 

e)  £LiZjiQ  (1200  ®C):  This  couple  exhibits  similar  behavior  to  the  Pi/Nb 
couple  analyzed  previously.  Here,  basically  the  penetration  depth  of 
the  elements  is  controlled  by  the  motion  of  Nb  in  the  Nb-rich  side, 
which  has  the  higher  solidus  temperature. 

0  WNb  (1400  °C);  The  behavior  of  this  couple  is  similar  to  that 
observed  in  the  same  system  at  1200  °C,  except  that,  as  expected, 
larger  penetration  depths  for  each  component  are  obtained, 
e)  (1400  ®C):  At  1400  °C  the  mobility  of  A1  is  higher,  and 

consequently  the  penetration  depth  of  it  is  larger.  At  this 
temperature  Ti  penetrates  deeper  into  Pi  side,  however,  its  draw 
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depth  from  the  P3  side  is  the  same  as  that  at  1200  ®C.  As  mentioned 
previously  the  draw  depth  of  Nb  from  the  P3  side  indeed  decreases 
with  increasing  temperature.  This  observatation  and  the  absence  of 
maxima  and  minima  in  the  composition  profile  of  Ti  at  1400  ®C  and 
its  almost  equal  depth  of  draw  with  that  of  A1  penetration  on  the  P3 
side,  contrary  to  the  pattern  at  1200  ®C,  suggest  that  at  the  latter 
diffusion  diffusion  temperature,  uphill  diffusion  of  Nb  have  occurred. 
Further  inve  igation  is  required  to  understand  the  composition 
patterns  in  this  system. 


The  composition  paths  for  the  different  diffusion  couples  are 
presented  in  3  parts  depending  on  their  position  over  the  Gibbs 
triangle: 

a)  Nb-rich  corner 

b)  Constant  Nb:Ti  ratio  zone 

c)  Ti-rich  corner 


5.4.1  The  Nb-Rich  Corner 


The  composition  paths  of  the  diffusion  couples  involving  a  pure 
Nb  side  (couples  #8  and  #6)  show  an  interesting  similarity  to  the  line 
joining  the  end  compositions  of  the  couple  as  shown  in  Figures  5.38 
and  5.39  for  interdiffusion  temperatures  of  1200  ®C  and  1400°C, 
respectively.  This  composition  path  behavior  indicates  [69]  equal 
mobility  of  the  elements  in  this  corner,  consistent  with  the  results 
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shown  in  Table  5.1.  The  composition  path  for  the  other  couple  in  this 
corner,  #11  (Figure  5.40),  shows  a  visible  bending  of  the  path  with  a 
crossover  point  close  to  the  ^2  composition.  In  this  couple  the  A1 
concentration  gradient  is  almost  zero,  however,  as  shown  in  Figure 
5.29  there  is  an  apparent  redistribution  of  A1  possibly  due  to 
unequal  mobility  of  Ti  and  Nb  in  this  couple. 

5.4.2  The  Constant  Nb:Ti  Ratio  Zone 

For  this  section  the  composition  paths  for  diffusion  couples  #10, 
#9  and  #7,  which  are  given  in  Figures  5.41,  5.42  and  5.43 
respectively,  are  considered.  The  similarity  of  the  composition  paths 
for  these  3  couples  is  evident.  Both  ends  of  the  composition  paths  are 
directed  toward  the  Ti  corner  and  the  midsection  at  the  crossover 
seems  to  follow  a  path  parallel  to  a  constant  Ti  composition. 

5.4.3  The  Ti-Rich  Comer 


There  is  only  one  composition  path  for  the  discussion  of  the 
diffusivity  behavior  in  this  corner  (couple  #12).  The  pattern  followed 
by  this  path  is  somehow  similar  to  that  obtained  for  couple  #  11,  but 
with  a  more  pronounced  tendency  to  deviate  from  the  line 
connecting  the  composition  ends  and  it  bends  toward  the  Nb-Ti  side 
of  the  couple  as  shown  in  Figure  5.44.  A  high  magnification  of  the 
path  near  the  Ti  comer  is  shown  in  Figure  5.45.  It  is  observed  from 
these  Figures  that  the  end  of  the  paths  at  both  ends  point  toward  the 
AI  comer  of  the  Gibbs  triangle,  which  implies  that  A1  moves  faster 
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Ti-rich  side  of  the  composition  path  for  diffusion 
couple  #12  (Pl/Ti). 
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than  Nb.  The  relative  mobility  of  Ti  and  A1  may  be  deduced  by  the 
displacement  of  the  crossover  point  toward  the  Ti  cornei  suggesting 
a  relative  higher  mobility  of  Ti.  This  matter  will  be  discussed  later  in 
detail  (section  5.6). 


Kirkendall  Shifts 


Two  different  methods  were  employed  to  measure  the 
Kirkendall  shift  in  the  systems  under  study.  The  first  of  them  was 
done  experimentally  by  measuring  the  distance  between  the  two 
interfaces  present  in  the  sandwich-like  couples  before  and  after  the 
diffusion  treatment.  The  resultant  difference  was  then  divided  by  2 
to  get  the  net  displacement  of  the  markers.  The  second  method 
employed  the  same  computer  program  used  for  the  evaluation  of  the 
Matano  Interface  in  the  binary  system.  In  this  case,  the  Matano 
plane  determination  was  conducted  over  the  profiles  of  the  3 
components  for  each  diffusion  couple.  In  general,  the  calculated 
Matano  interface  position  was  approximately  the  same  for  the  three 
elements.  Table  5.2  shows  the  Kirkendall  shift  as  calculated  by  the 
above  mentioned  methods.  There  is  an  excellent  agreement  between 
the  two  methods,  except  for  the  couple  #11  (P2/Po)>  for  which  both 
methods  show  similar  displacements  but  in  opposite  directions.  It 
was  observed  that  the  agreement  between  results  of  both  methods  is 
better  when  the  computer  program  was  applied  to  the  experimental 
A1  composition  profiles.  The  reason  for  this  could  be  the  absence  of 
critical  points  in  the  A1  profiles  in  most  cases.  One  problem  associated 
with  the  experimental  method  was  the  lack  of  statistical 
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measurements  as  only  a  few  number  of  yttria  markers  were 
available.  For  the  case  of  couple  #12,  they  were  not  even  found  as 
they  were  probably  dissolved  in  the  Ti  matrix.  Instead,  in  this  case 
the  measurements  were  made  using  oxide  markers  developed  in  the 
bonded  interface  during  the  hot-pressing  operation. 

At  1200  ®C,  the  diffusion  couple  with  pure  Ti  (#12)  shows  the 
largest  displacement  observed  with  a  Kirkendall  shift  toward  the  Ti 
side  of  the  couple.  Regarding  the  two  1400  ®C  treated  couples,  a 
significant  marker  displacement  (Kirkendall  shift)  was  observed  for 
couple  #6  pointing  toward  the  ternary  side  of  the  couple  (rich  Ti, 
poor  Nb).  As  to  the  other  couple  treated  at  1400  ®C  (#7)  and  those 
interdiffused  at  1200  ®C,  the  Kirkendall  shift  is  small  and  in  some 
cases  negligible.  For  example,  the  P1/P3  couples  with  constant  Nb:Ti 
ratio  and  15  at%  A1  (couples  #7  and  #9)  exhibited  a  negligible 
displacement  of  the  markers  at  either  temperature,  1400  °C  or  1200 
®C.  Assuming  only  the  A1  composition  gradient,  then  one  would 
expect  a  larger  shift  in  the  P1/P3.  couple  than  in  the  I33/P2  couple 
(5at%  AI).  However,  the  consideration  of  the  increase  in  Nb  and  Ti 
gradients  suggests  that  more  Nb  and  Ti  flow  into  the  ternary  side  in 
the  P3/P1  couple  than  they  do  in  the  P3/P2  couple. 

5.6  Crossover  Displacements 

As  mentioned  in  Chapter  2,  the  crossover  point  of  a 
composition  path  in  the  Gibbs  triangle  is  displaced  from  the  mid 
point  of  the  line  connecting  the  end  compositions  when  the 
penetration  tendencies  of  the  components  are  concentration- 


dependent.  The  composition  of  the  crossover  point  should  correspond 
to  the  composition  of  the  original  weld  interface,  which  is  given  by 
the  composition  of  the  Matano  interface.  In  Figures  5.38  to  5.44  the 
Matano  interface  compositions  are  marked  on  the  composition  paths. 
There  is  a  good  correlation  between  the  calculated  and  crossover 
compositions,  mainly  in  the  couples  involving  Pi,  P2.  and  P3  alloys. 

5.6  Discussion 


The  discussion  of  the  diffusivity  in  the  ternary  Nb-Ti-Al  p-solid 
solution  is  based  on  the  results  obtained  for  penetration  tendencies 
(see  Table  5.1),  composition  paths  (see  Figures  5.38-5.45),  and 
Kirkendall  shifts  (see  Table  5.2).  An  analysis  of  the  composition 
profiles  permitted  to  evaluate  the  approximate  penetration  ratios  of 
the  components  in  the  different  couples.  These  ratios  are  given  in 
Table  5.3,  and  it  can  be  observed  that  they  depend  strongly  on  the 
composition  of  the  alloys.  The  discussion  of  these  results  will  be 
focused  on  three  different  regions  and  emphasis  will  be  placed  on 
the  diffusion  composition  patterns  following  DeHoff  et.  al.  [69]. 
a)  Nb/Ti  Constant  Ratio 

The  composition  paths  of  the  3  couples  in  this  region  (#10,  #9, 
#7)  as  shown  in  Figures  5.41  to  5.43  are  similar  in  behavior.  It  is 
observed  that  the  end  segments  are  directed  toward  the  Ti  corner, 
therefore,  indicating  that  Ti  is  the  fast  diffuser.  The  midsection  of 
these  paths  are  parallel  to  the  Nb-Al  side  suggesting  that  Al  and  Nb 
have  similar  penetration  tendencies.  When  the  midsections  of  the 
P3/P1  couples  annealed  at  two  temperatures  are  compared  in  more 
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detail  (#7  and  #9),  it  is  observed  that  the  midsection  of  couple  #7 
(1400  ®C)  presents  a  slight  rotation  toward  the  A1  corner,  suggesting 
that  the  mobility  of  Nb  increases  as  the  temperature  increases.  It  is 
also  observed  that  deviation  of  the  composition  path  from  the 
straight  line  behavior  in  couples  #7  and  #9  decreases  with  increasing 
the  annealing  temperature,  which  suggests  that  the  mobility  of  A1 
and  Nb  is  closer  to  that  of  Ti  at  higher  temperatures. 

As  far  as  the  penetration  tendency  ratios  are  concerned  (Table 
S.3),  the  three  couples  show  the  highest  penetration  ratio  for  the 
element  Ti,  consistent  with  the  observations  regarding  the 
composition  paths.  However,  when  the  ratios  for  Nb  and  A1  are 
compared,  an  artificially  large  penetration  of  Nb  is  found  for  the 
P1/P3  couple  (#9)  annealed  at  1200  ®C.  As  demonstrated  by  DeHoff  et. 
al.  [69],  the  penetration  depths  calculated  from  the  microprobe 
analysis,  particularly  for  elements  that  are  not  indeed  mobile,  do  not 
represent  the  physical  distance  over  which  the  elements  have 
moved.  In  this  case  the  imbalance  in  distribution  of  Al  and  Ti  alters 
the  percentage  of  Nb,  This  is  supported  by  the  facts  that  Nb  shows  a 
larger  depth  at  1200  °C  than  it  does  at  1400  ®C  in  the  P3  side,  and 
that  the  reaction  depth  of  Ti  at  1200  °C  is  much  larger  than  that  of 
Al. 

A  comparison  of  the  penetration  depths  for  P1/P3  and  P2/P3 
couples  (Table  5.1)  reveals  that  in  general  the  depths  are  much 
larger  for  the  P1/P3  couple,  which  had  a  larger  concentration  gradient 
of  the  elements.  The  penetration  depth  of  Al  on  the  P3  side  (Nb-50Ti) 
decreases  from  72  pm  to  31  pm  with  a  decrease  in  concentration 
gradient.  It  is  interesting  to  note  that  the  Nb  and  Ti  are  apparently 


drawn  from  much  deeper  in  the  ^1/33  couple,  however,  in  the  ^2/^3 
couple  their  depth  01  draw  is  similar  to  the  penetration  depth  of  Al. 
The  larger  depth  of  draw  of  Ti  in  the  P1/P3  couple  can  be  attributed 
to  the  larger  concentration  gradient  of  this  element  in  comparison  to 
the  |32/p3  couple.  Consistently,  in  both  couples  Nb  penetrates  to  a 
shorter  depth  into  the  counterpart  side  of  the  P3  than  does  Ti. 

As  far  as  Kirkendall  shift  for  these  couples  are  concerned,  the 
displacements  were  measured  to  be  very  small  but  always  toward 
the  ternary  side,  indicating  that  physically  more  Al  atoms  were 
transferred  to  the  P3  side.  This  observation  is  consistent  with  the 

larger  concentration  gradient  of  Al  than  those  of  Ti  and  Nb  for  these 

couples. 

b)  TI  Corner 

In  consistency  with  the  assumption  that  Ti  moves  much  faster 
than  Al  and  Nb  the  composition  path  of  p3/Ti  couple  (couple  #12) 
showed  the  tendency  toward  a  straight  line  behavior  (see  Figure 
5.44).  The  end  segment  at  the  P3  composition  is  directed  toward  the 
Al  corner  indicating  that  in  this  corner  Al  moves  faster  than  Nb. 
Consequently  the  composition  path  deviates  slightly  from  the 
straight  line  behavior,  but  due  to  a  much  faster  diffusivity  of  Ti  it 
crosses  over  very  close  to  the  Ti  corner. 

The  system,  in  general,  shows  the  highest  penetration  depths  in 
this  corner  consistent  with  the  fact  that  Ti  has  the  lowest  melting 
point  among  all  constituent  alloys  used  in  this  study.  As  shown  in 
able  5.1  all  elements  show  a  larger  depth  of  diffusivity  in  the  Ti 
side  than  in  the  P  j  side  consistent  with  the  high  solidus  temperature 
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of  the  latter  alloy.  Nb  tends  to  show  shorter  draw  and  penetration 
depths  than  Ti  and  A1  do. 

The  Kirkendall  Shift  measured  in  this  couple  is  also  one  of  the 
largest  observed  in  all  couples,  and  accordingly,  it  predicts  a  larger 
flux  of  Ti  into  the  pi  side  of  the  couple, 
c)  Nb  Corner 

The  penetration  ratio  values  (1:1:1)  given  in  Table  5.3  for 
couples  #6  and  #8  are  in  complete  agreement  with  the  matching  of 
composition  paths  to  the  straight  line  connecting  the  end 
compositions  observed  for  these  couples.  This  behavior  can  be 
interpreted  as  a  result  of  the  low  jump  frequencies  of  all  three 
components  in  this  region  of  the  diagram,  which  again  can  be 
correlated  to  the  solidus  temperature  of  the  alloys  involved  and  is 
supported  by  the  depth  of  penetration  measured  on  the  two  sides  of 
the  couples.  These  penetrations  consistently  show  a  very  small 
penetration  into  the  Nb  side  and  it  is  the  same  for  the  three 
elements.  The  penetration  depth  in  the  pi  side,  on  the  other  hand 
shows  a  deeper  diffusion  region  due  to  the  lower  solidus 
temperature  of  this  alloy.  The  low  value  of  the  Kirkendall  shift 
observed  at  1200  ®C  also  shows  that  the  mobility  of  the  three 
elements  is  similar.  However,  when  the  temperature  is  increased  to 
1400  ®C,  the  Kirkendall  shift  measured  is  very  large  (see  Table  5.2), 
which  is  difficult  to  interpret  since  it  indicates  that  a  large 
displacement  of  Ti  and  A1  atoms  toward  the  Nb  had  occurred.  Such 
observations  do  not  agree  with  the  composition  path  and  the 
penetration  depth  results.  Further  investigation  is  required  to 
resolve  these  inconsistencies. 
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The  results  of  this  study  suggest  that  the  diffusivity  of  the 

components  in  general  is  larger  when  the  amount  of  titanium  present 
in  a  couple  increases.  The  Kirkendall  shift  and  crossover  point 
displacement  results  also  indicate  that  the  mobility  of  the  system 

increases  as  the  titanium  content  increases,  but  they  also  suggest 
that  the  aluminum  owns  a  larger  mobility  than  niobium  as  both  the 
Kirkendall  shift  and  crossover  displacement  point  opposite  to  the 
niobium  comer.  In  accordance  with  these  results,  titanium  owns  the 
largest  jump  frequency  followed  by  aluminum  and  the  niobium. 

Another  observation  of  this  study  is  that  the  Kirkendall  and 

crossover  shifts  are  in  agreement  as  to  the  direction  of  the 

displacement,  and  the  composition  at  the  Matano  interface  correlates 
well  with  the  crossover  composition.  The  mobility  of  all  components 
is  larger  close  to  the  Ti  corner  of  the  diagram.  The  jump  frequency  of 
the  system  changes  more  drastically  as  going  from  the  Nb  comer  side 
of  the  system  toward  the  Ti  comer.  This  behavior  causes,  therefore,  a 
similar  crossover  shift  in  the  couples  involving  end  compositions  in 
the  Ti-rich  section  of  the  ternary  system,  and  a  larger  crossover  shift 
in  the  couples  involving  Nb-rich  end  compositions. 

5.7  Summary 

Contrary  to  the  behavior  of  Nb-4.5  at%  A1  solid  solution,  no 
propensity  to  internal  oxidation  was  detected  in  the  ternary  alloy. 
This  observation  suggests  that  the  addition  of  Ti  reduces  the 
susceptibility  of  the  alloys  to  internal  formation  of  AI2O3.  Despite  the 
observation  of  two-phase  microstructures  at  room  temperature,  it 
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CHAPTER  6 
CONCLUSIONS 

The  diffusion  behavior  in  the  binary  Nb-Al  and  the  ternary  Nb- 
Ti-Al  systems  were  studied  employing  Matano-Boltzmann  diffusion 
couple  method.  The  results  obtained  have  led  to  the  following 
conclusions: 

1-  The  binary  Nb-4.5at%  A1  alloy  was  very  susceptible  to  internal 
oxidation.  However,  the  addition  of  Ti  to  the  b.c.c.  solid  solution 
completely  retarded  the  internal  oxidation. 

2-  The  zirconia  particles  used  as  markers,  were  unstable  in  the  Nb-Al 
p  solid  solution  and  iecomposed  at  high  temperatures.  The 
substitution  of  the  markers  with  yttria  proved  that  this  oxide  is 
completely  inert  in  the  binary  Nb-Al  ^-solid  solution,  however,  in  the 
presence  of  Ti  some  reactivity  was  observed. 

3-  The  interdiffusion  coefficient  of  the  binary  Nb-Al  system  varied 
linearly  with  composition  in  the  range  of  1-3  at%  A1  and  it  showed  a 
linear  Arrhenius  behavior  as  a  function  of  temperature.  The  average 
interdiffusion  coefficient  in  the  temperature  range  of  1350  ®C  to 
1550  °C  can  be  expressed  through  the  following  equation: 

D  =  0.1136  exp  (-80,105  (cal/mole)/RT)  cm^/sec 
The  activation  energy  was  independent  of  the  composition  in  the  1-3 
at%  range.  The  average  intrinsic  diffusivity  of  A1  in  the  same 
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temperature  range,  calculated  based  on  a  dilute  solution  assumption 
is  given  by: 

Dai  =  0.1161.  exp  (-80,165  (cal/mole)/RT)  cm^/sec 
4-  Kirkendall  shift  measurements  correlate  well  with  the  calculated 
crossover  displacements  in  the  ternary  Nb-Ti-Al  system.  The 
direction  and  magnitude  of  displacements  in  conjunction  with  an 
analysis  of  the  composition  paths  suggest  that  Ti  is  the  fastest 
component  of  the  system  at  all  compositions,  with  a  very  small 
difference  in  the  relative  mobilities  of  A1  and  Nb.  The  addition  of  Ti 
increases  the  general  diffusivity  of  the  system  partially  due  to  a 
reduction  in  the  solidus  temperature  of  alloys. 

The  occurrence  of  maxima  and  minima  in  the  composition  profiles 
obtained  by  microprobe  analysis  correlates  with  the  existence  of 
small  composition  gradients  or  very  low  relative  mobilities. 

6-  Results  of  this  study  show  that  the  3  elements  Nb,  Ti  and  Al,  have 
similar  but  very  low  jump  frequencies  in  the  Nb  -rich  corner  of  the 
Gibbs  triangle.  In  the  Ti-rich  corner  the  jump  frequencies  of  these 
elements  are  unequal  but  much  higher  than  the  Nb  jump  corner.  Ti 
has  the  highest  jump  frequency  in  this  comer. 
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APPENDIX  A 

PROBLEMS  ASSOCIATED  WITH  THE  MEASUREMENT  OF  COMPOSITION 
PROFILES  BY  MICROPROBE  ANALYSIS 

The  evolvement  of  electronic  analysis  techniques  such  as  the 
Electron  Microprobe  for  the  study  of  diffusivity  process  have 
provided  the  capability  to  measure  gradient  concentrations  at  very 
low  spatial  resolutions  (1pm)  that  traditional  techniques  would  have 
never  been  able  to  reach.  This  has  permitted  to  study  diffusion  at 
very  short  scales,  such  as  low  temperature  and/or  short  time 
treatments.  However,  the  success  of  the  application  of  this  technique 
relies  on  the  accuracy  of  the  data  obtained.  In  order  to  report 
reliable  results,  several  composition  profiles  were  measured  in  each 
of  the  diffusion  couples  studied  in  this  research.  At  least  two 
matching  profiles  were  obtained  for  each  couple  (Figures  A.l),  and 
the  one  with  the  smoother  curve  was  chosen  for  calculating  of  the 
Matano  plane  position. 

The  inconsistency  in  composition  profiles  obtained  for  an 
element  in  the  same  couple  was  found  to  occur  in  two  ways. 
Sometimes  a  constant  shift  in  the  composition  profile  was  detected. 
Due  to  a  change  in  parameters  associated  with  the  calibration  of  the 
instrument,  the  entire  profile  was  shifted  with  respect  to  the  initial 
profile.  Since  this  displacement  was  equal  for  all  compositions,  no 
effect  on  the  calculated  depth  of  penetration  was  '  detected.  The 
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Composition  profiles  of  A1  after  interdiffusion 
treatment  of  couple  #9  at  1200  °C  for  6  hours. 
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computation  of  the  Matano  interface  was  not  affected  either,  since 
the  end  compositions  employed  in  such  calculations  were  those 
obtained  while  measuring  the  resulting  composition  profiles.  Figure 
A2  shows  this  behavior  for  the  Nb  composition  profiles  in  couple  #9. 
The  other  effect  on  the  measurement  of  the  composition  profiles  was 
that  due  to  the  existence  of  markers  close  to  the  point  from  which 
the  profile  was  obtained.  Figure  A3  shows  the  BSE  micrograph  of  a 
diffusion  profile  measured  in  diffusion  couple  #10  and  the  presence 
of  yttria  markers.  When  the  microprobe  analysis  of  the  couple  was 
performed,  the  closeness  of  the  markers  produced  compositions 
profiles  as  that  shown  in  Figure  A4.a,  When  the  analysis  was 
performed  in  zones  far  enough  from  the  markers,  normal  profiles  as 
that  shown  for  the  same  couple  in  Figure  A4.b  was  obtained. 
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Composition  profile  of  Nb  after  interdiffusion 
treatment  of  couple  #9  at  1200  ®C  for  6  hours. 
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_  10  urn 

Figure  A3  BSE  micrograph  showing  microprobe  analysis 

contamination  close  to  yttria  markers  in  couple 

#10  (P3/P2). 
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Figure  A4  Composition  profiles  of  Nb  after  interdiffusion 

treatment  of  couple  #10  at  1200  ®C  for  6  hours. 

a)  Measured  close  to  yttria  markers 

b)  Measured  away  from  markers. 
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